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The  feasibility  of  developing  a  diode  laser  Doppler  velocimeter  (DLDV)  for 
oceanographic  measurements  of  current  and  turbulence  has  been  demonstrated  by 
Flow  Industries,  Inc.  (FLOW)  in  a  DOD  SBIR  Phase  I  investigation  sponsored  by 
the  Office  of  Naval  Research  (see  Phase  I  final  report,  Schedvin  and  Liu,  1984 
—  referred  to  as  [I]  hereafter).  The  Phase  II  work,  which  is  documented  in 
this  report,  was  directed  toward  the  development  of  a  complete  DLDV  sensor 
system  that  could  be  easily  adapted  to  a  number  of  field  measurement  applica¬ 
tions.  This  required  the  development  of  both  single-component  and  multi- 
component  DLDV  sensors  and  the  development  of  a  flexible  data  analysis 
electronics  system  that  can  be  used  under  a  variety  of  flow  conditions.  A 
direct  product  of  the  Phase  II  effort  is  an  oceanic  boundary  layer  measurement 
system  using  high-resolution  sensor  clusters  (HRCs)  .  The  DLDV  sensors 
developed  for  the  measurement  system  serve  as  prototypes  of  the  commercial 
sensors  to  be  manufactured  and  marketed  in  Phase  III. 

Part  of  the  Phase  II  work  also  involved  supporting  the  development  and 
field  deployment  of  moderate-resolution  clusters  (MRCs),  each  consisting  of 
three  small,  partially  ducted  current  meters  and  a  fast-response  temperature 
and  conductivity  sensor  (McPhee  and  Smith,  1976).  The  spatial  resolution  of 
the  MRC  is  limited  to  about  10  to  30  cm  by  the  response  of  the  current  meters 
and  colocation  of  the  sensors.  A  detailed  description  of  the  MRC,  the  field 
deployment  and  results  is  given  elsewhere  (McPhee,  1984,  1985a). 

The  Phase  II  effort  focused  mainly  on  the  further  development  of  the  DLDV 
sensor.  The  Phase  I  results  (l]  indicated  the  feasibility  of  an  underwater 
DLDV  sensor.  However,  a  number  of  developmental  steps  were  still  needed  to 
produce  an  instrument  of  commercial  interest  for  oceanographic  measurements. 
The  two  primary  objectives  of  the  Phase  II  work  are  described  below. 


(1)  Development  of  a  compact,  reliable  DLDV  sensor  with  optional  flow 
direction  sensing. 

This  objective  represented  the  main  developmental  effort  of  Phase  II. 
Primary  tasks  were  the  development  of  Doppler  signal  analysis  electronics 
and  the  selection  and  implementation  of  the  optimal  method  of  flow 
direction  sensing.  Other  tasks  included  the  optimization  of  the  laser 
diode  optical  configuration,  the  photodetector,  scattered  light  pickup, 
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and  undervater  packaging  for  low  flow  disturbance  and  ease  of  maintenance. 
The  result  is  a  high-resolution,  low-power,  compact  velocity  sensor  with 
many  oceanographic  applications. 

(2)  Development  of  a  high-resol  tion  turbulent  flux  measurement  cluster 
for  oceanic  boundary  layer  measurements . 

The  high-resolution  clusters  use  the  DLDV  sensors  to  provide  a  three- 
component  velocity  measurement  and  include  fast-response  temperature  and 
conductivity  sensors.  A  modular  approach  utilizing  the  single-component 
DLDV  sensor  and  associated  analysis  electronics  was  adopted.  A  standard 
digital  interface  bus  is  used  for  all  instruments,  with  microprocessor 
control  of  sampling  and  data  storage.  This  allows  great  flexibility  in 
the  auxiliary  sensor  and  the  configuration  of  the  DLDV  sensors. 

Our  goal  was  to  achieve  a  spatial  resolution  of  about  2  cm  for  the  HRCs, 

10  times  greater  than  that  of  the  MRCs .  This  represents  a  significant  advance 
in  oceanographic  instrumentation  for  boundary  layer  measurements  toward 
complete  resolution  of  the  turbulent  fluxes.  While  field  sensors  with  such  a 
high  resolution  are  lacking,  the  turbulent  fluxes  of  momentum,  heat  and  mass 
at  this  fine-scale  (and  high-frequency)  range  may  prove  to  contribute 
significantly  to  the  total  fluxes  in  the  oceanic  boundary  layer  under  ice, 
such  as  has  been  observed  in  the  atmospheric  boundary  layer  (Lumley  and 
Panofsky,  1964). 

The  modular  instrument  design  with  flexible  digital  interfacing  and 
control  will  provide  a  versatile  oceanographic  sensor.  While  the  initial 
commercial  prototypes  are  intended  for  oceanic  boundary  layer  measurements, 
the  instruments  could  also  be  used  on  towed  or  self-propelled  platforms  and  on 
vertically  profiling  packages.  ' 

Field  testing  was  an  important  part  of  the  Phase  II  effort  and  had  two 
primary  goals.  First,  the  field  tests  provided  an  assessment  of  the  relia¬ 
bility  of  the  DLDV  sensors  and  allowed  evaluation  of  sensor  performance  with 
oceanic  scatterer  concentrations.  This  information  was  used  to  improve  the 
design  of  the  sensor  and  the  Doppler  signal  analysis  electronics.  Second,  the 
field  tests  provided  the  environment  for  examining  the  capabilities  of  the 
moderate-  and  high-resolution  clusters  for  conducting  turbulence  flux 
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measurements  in  the  oceanic  boundary  layer  under  ice,  as  well  as  providing 
important  data  sets  for  Arctic  Ocean  boundary  layer  studies. 

This  report  documents  the  research  and  development  conducted  in  the 
Phase  II  investigation  that  led  to  the  completion  of  two  experimental  DLDV 
underwater  units  deployed  during  MIZEX  84  and  subsequently  the  completion  of 
three  prototype  velocity-temperature-conductivity  (V-T-C)  HRCs  deployed  during 
AIWEX  85.  Documentation  of  the  two  field  experiments,  MIZEX  84  and  AIWEX  85, 
and  some  of  the  field  results  are  also  included.  Section  2  describes  the 
individual  instruments  and  equipment  that  make  up  the  HRC.  In  Section  3,  we 
present  the  performance  characteristics  of  the  HRC  sensors.  The  software  for 
data  acquisition  and  analysis  is  given  in  Section  4.  Section  5  documents  the 
field  experiments  and  the  results  analyzed  to  date.  Finally,  a  summary  to 
assess  the  performance  of  the  HRC  in  the  Arctic  environment  and  recommendations 
for  future  work  are  presented  in  Sections  6  and  7,  respectively. 

2.  IMSTROMEHTATION  OF  THE  HIGH-RESOLOTION  CLUSTER 

The  high-resolution  cluster  consists  of  three  DLDV  underwater  units,  a 
high-frequency  thermistor  probe  and  a  dual-electrode  conductivity  probe 
mounted  on  an  octagonal  holder  made  of  stainless-steel  tubes  and  parts.  Most 
of  the  research  and  development  concentrated  on  the  DLDV  units.  During  the 
Phase  I  feasibility  demonstration,  we  constructed  and  tested  a  bench  test  model 
of  the  DLDV  [I]  .  Then,  in  Phase  II,  we  built  two  experimental  DLDV  units  and 
deployed  them  during  MIZEX  84  to  demonstrate  their  field  applicability.  Both 
a  counter  and  a  frequency  tracker  were  used  as  the  Doppler  signal  processors 
in  an  attempt  to  determine  the  optimal  processor  for  our  intended  field 
measurements.  Based  on  the  findings  obtained  in  the  laboratory  tests  and 
MIZEX  84,  we  redesigned  and  assembled  nine  prototype  DLDVs  and  integrated  them 
into  three  HRCs,  which  were  deployed  subsequently  during  AIWEX  85.  Figures  1 
and  2  show  a  schematic  and  photographs  of  the  HRC. 

The  thermistor  probes  and  conductivity  probes  selected  for  the  HRCs  have 
been  used  in  oceanographic  measurements.  Their  performance  characteristics 
have  been  reported  in  the  literature  (Meagher  et  al.,  1982).  Several  new 
features  were  incorporated  into  the  circuit  design  to  be  compatible  with  our 
measurement  strategy. 
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Schctnatic  of  the  HRC  Consisting  of  a  Three-Axis  DLHV  Syste 
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a.  Basic  Building  Block  of  the  Transmitting  Module 


b.  Layout  of  the  HRC  Sensors 


Figure  2,  Photographs  of  the  High-Resolution  Cluster 
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The  sensors  on  each  HRC  were  connected  to  a  processor  developed  and 
assembled  by  FLOW,  which  provides  power  to  drive  the  sensors,  provides  signal 
amplification,  filtering  and  preprocessing,  and  interfaces  the  signal  outputs 
with  the  host  computer. 

In  this  section,  we  describe  in  detail  the  individual  sensors.  First  we 
discuss  the  DLDV  modules,  then  the  temperature/conductivity  probe  pair,  and 
finally  the  HRC  processors. 

2.1  DLDV  Modules 

The  concept  of  the  DLDV  and  the  design  considerations  are  given  in  detail 
in  [l]  .  In  this  report,  we  describe  the  actual  design  and  assembly  of  the 
prototype  units.  The  DLDV  consists  of  a  transmitting  module  and  a  receiving 
module.  A  schematic  of  the  two  modules  set  up  for  velocity  measurements  to¬ 
gether  with  detailed  drawings  is  shown  in  Figure  3  (also  see  Figure  14  of  [I]). 

2.1.1  Transmitting  Module 

As  shown  in  Figure  3,  the  transmitting  module  includes  the  electro-optical 
(EO)  components  to  form  a  diode  laser,  to  split  the  laser  beam  and  to  configure 
Che  beams  into  the  fringe  mode  of  arrangement  for  use  with  laser  velocimetry. 

The  EO  components  are  mounted  on  an  aluminum  chassis,  which  is  in  turn  placed 
in  an  underwater  housing  made  of  a  stainless-steel  tube  (15.2  cm  long  and 
3.8  cm  OD) .  The  diode  laser  consists  of  a  5-mW  laser  diode  (Hitachi  HL7801G) 
mounted  on  a  three-dimensional  miniature  stage,  a  collimator  (Fischer  40x  objec¬ 
tive),  a  beam  splitter  and  a  focal  lens  (8.5  focal  length  in  air).  The  nominal 
wavelength  of  the  diode  laser  is  793  nm,  which  is  in  the  deep  red  regime  fnot 
visible  unless  viewed  with  a  fluorescent  screen  (Optical  Engineering)]  . 

The  beam  splitter  is  formed  by  gluing  five  1-cm  prisms  manufactured  by 
Melles  Griot.  Two  of  the  prisma  are  preassembled  by  the  manufacturer  (Melles 
Griot)  to  form  a  cube  beam  splitter  coated  for  laser  diode  applications.  The 
rest  of  the  prisms  are  glued  to  the  cube  beam  splitter  using  an  optical 
adhesive  that  is  cured  with  an  UV  light.  A  special  jig  mounted  on  an  optical 
table  is  used  to  align  the  prisms,  and  the  alignment  is  monitored  using  a  5-mW 
He-Ne  laser  during  the  entire  curing  period.  The  error  in  parallelism  of  the 
two  split  beams  is  maintained  to  be  within  ^  5  millidegrees  in  both  planes. 

The  prisms  and  the  focal  lens  are  coated  with  a  HEBBAR  anti-reflection  coating 
(Melles  Griot)  to  reduce  reflection  and  to  enhance  optical  transmission. 
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Figure  3.  Design  of  the  DLDV  Transmitting/Receiving  Modules 


To  maintain  a  constant  laser  power  output  and  to  protect  the  laser  diode 
from  large  electrical  surges  or  spikes  present  in  the  power  source,  a  special 
circuit  (Figure  4)  was  designed  and  inserted  between  the  diode  and  the  DC 
power  supply.  The  desired  forward  current  may  be  set  using  a  10-turn  pot.  A 
feedback  loop  using  the  output  from  a  built-in  photodiode,  which  monitors  the 
laser  power  output,  controls  the  forward  current  to  maintain  a  constant  laser 
power  output.  For  example,  when  the  surrounding  temperature  drops,  the  diode 
laser  becomes  more  efficient.  The  feedback  loop  automatically  reduces  the 
forward  current  to  maintain  a  constant  laser  power. 


Figure  4.  Feedback  Circuitry  for  Control  of  the  Laser  Power  Output 


2.1.2  Receiving  Module 

The  receiving  module  consists  of  the  EO  components  for  collecting  and 
sensing  the  scattered  light  from  tracer  particles  passing  through  the  focal 
volume  of  the  laser  beams.  A  detailed  schematic  of  the  module  is  shown  in 
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Figure  3.  The  module  is  enclosed  in  an  underwater  housing  made  of  a  stainless- 
steel  tube  (10.2  cm  long  and  3.8  cm  OD) .  A  3-mm-diameter  fused  fiber  optic 
rod  is  used  as  the  collector  of  scattered  light.  The  advantage  of  using  a 
fiber  optic  rod  is  twofold:  it  avoids  the  necessity  for  precision  alignment, 
and  it  reduces  light  attenuation  by  bringing  the  rod  end  reasonably  close  (2 
to  3  cm)  to  the  focal  volume  without  introducing  significant  flow  disturbance. 
Scattered  light  falling  on  the  rod  end  is  internally  reflected  by  the 
individual  fibers  (about  70,000  in  the  rod  used)  and  routed  to  the  photodiode 
(Devar  Inc.,  Type  539-01-5,  with  a  5-mm  sensing  area).  The  driver  of  the 
photodiode  is  mounted  on  a  printed  circuit  board  next  to  the  diode.  In 
addition  to  driving  the  photodiode,  the  circuit  provides  the  functions  of 
low-pass  filtering  and  amplification  (up  to  a  gain  of  20)  of  the  signals.  The 
output  signals  and  the  DC  power  (+  15  volts)  are  transmitted  to  and  from  the 
surface  through  the  main  underwater  cable  with  a  pigtail  (branched  out  from  a 
junction  box)  plugged  into  a  five-pin  connector  at  the  end  of  the  module. 

2 .2  Temperature/Conductivity  Probe  Pair 

Figure  5  shows  a  drawing  of  the  temperature/conductivity  probe  pair  used 
on  the  HRC  cluster.  It  consists  of  a  thermistor  probe  (FASTIP,  Model  FP07 , 
Thermometries,  Inc.)  and  a  dual-electrode  conductivity  probe  glued  together 
with  the  three  tips  3  mm  apart  and  forming  a  triangle.  The  thermistor  is 
housed  at  the  tip  of  a  shock-resistant  glass  tube.  The  small  bead  thermistor 
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Figure  5.  Schematic  of  Temperature/Conductivity  Probe  Pair 
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has  an  extremely  thin  glass  coating.  Its  response  time  is  7  ms  in  water  as 
measured  in  a  plunging  test.  The  nominal  resistance  of  the  thermistors  is 
10  kohm  at  25®C.  The  conductivity  probe  consists  of  two  thin  platinum  wires 
(0.05  cm  in  diameter)  each  sealed  in  a  thin  glass  tube  with  the  tip  exposed 
and  platinized  with  a  platinum  black  solution  manufactured  by  Yellow  Spring 
Instrument  Co.  The  probe  is  similar  to  that  reported  by  Meagher  et  al.  (1982) 
with  the  wires  2.5  times  smaller  in  diameter  than  their  probes. 

The  selection  of  the  dual-electrode  probe  rather  than  the  larger  conduc¬ 
tivity  cells,  such  as  those  of  the  Neil  Brown  or  Sea-Bird  cells,  is  due 
primarily  to  the  compatibility  of  the  DLDV  sensors  in  terms  of  physical  size, 
configuration,  and  dynamic  response.  Our  goal  is  to  achieve  an  overall 
spatial  resolution  of  1  to  2  cm  for  the  V-T-C  high-resolution  cluster. 
According  to  Meagher  et  al .  (1982),  the  dynamic  response  of  a  dual-electrode 
probe  (0.05-inch-diameter  platinum  wires)  is  about  100  cpm  (cycles/meter)  at 
the  3-4^  point.  Our  probe  should  have  better  response  because  the  size  of  the 
wire  is  about  2.5  times  smaller  than  theirs.  The  trade-off  in  using  the 
dual-electrode  probe  is  some  sacrifice  of  calibration  stability.  Meagher 
et  al.  (1982)  measured  a  drift  rate  of  this  type  of  probe  to  be  approximately 
0.01  S/m  (siemens/meter;  1  S/m  =  0.1  {jrfio/cm) .  Such  a  slow  drift  has  little 
effect  on  the  density  flux  measurements  because  the  differential  rather  than 
the  absolute  conductivity  is  of  interest  to  us. 

It  is  interesting  to  point  out  that  the  dual-electrode  probes  are  less 
prone  to  noncatastrophic  fouling  than  the  large,  more  stable  conductivity 
cells.  Experience  with  three  types  of  probes  towed  in  the  seasonal  thermocline 
of  the  Sargasso  Sea  shows  that  the  1.5-cm  Neil  Brown  cell,  the  Neil  Brown 
planar  cell,  and  the  dual-electrode  probes  were  fouled  35%,  42%,  and  8%  of  the 
time,  respectively  (Meagher  et  al.,  1982). 

2 .3  HRC  Proceasore 

The  function  of  the  processor  is  to  provide  DC  and  AC  power  to  drive  the 
various  electronic  active  elements  of  the  DLDV  and  temperature/conductivity 
sensors  and  to  receive  and  condition  (filtering,  amplification  and  DC  voltage 
offsetting)  the  raw  signals  output  from  these  sensors.  The  conditioned,  analog 
data  are  then  converted  into  digital  form  to  be  recorded  with  a  COMPAQ  computer 
using  a  12-bit  D/A  converter.  Appendix  A  is  a  set  of  operating  instructions 
for  the  software  controlling  the  HRC  processors. 
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2.3.1  DLDV  Processor 

An  important  part  of  the  Phase  II  effort  was  the  development  of  flexible 
data  analysis  electronics  for  the  Doppler  signal.  Standard  approaches  to 
Doppler  signal  analysis  are  the  use  of  period  counters  and  frequency  trackers. 
Because  the  output  power  of  laser  diodes  is  fairly  low,  small  sampling  volumes 
are  necessary  for  adequate  scattered  light  intensity.  This  requirement, 
combined  with  the  low  to  moderate  scatterer  densities  expected  in  the  open 
ocean  and  in  particular  in  the  Arctic,  indicates  that  a  counter  processor  is 
probably  the  best  choice  for  the  Doppler  signal  analysis  electronics  (Krumholz 
and  Murphy,  1974).  During  MIZEX  84,  we  used  both  a  tracker  and  a  counter  to 
process  the  Doppler  signals.  The  superiority  of  the  counter  to  the  tracker 
was  confirmed  from  the  MIZEX  84  results  (see  Section  5.1).  Ther  fore,  we 
decided  to  adopt  the  counter  approach. 

The  approach  used  in  a  counter  processor  is  to  first  band-pass  filter  the 
photodetector  output  to  remove  the  pedestal  signal.  This  is  a  low-frequency 
component  that  is  associated  with  the  variation  in  the  average  scattered  light 
intensity  as  a  particle  moves  through  the  sampling  volume.  The  remaining 
Doppler  frequency  burst  (Figure  6)  is  then  fed  to  a  period  counter  processor. 
This  processor  determines  the  time  interval  for  a  number  of  zero  crossings  of 
the  Doppler  signal.  The  time  interval  is  obtained  by  counting  a  high-frequency 
clock.  Also,  valid  zero  crossings  are  determined  by  simultaneously  checking  to 
see  that  the  Doppler  signal  crosses  an  appropriate  threshold  (either  above  or 
below  0  volts)  after  each  zero  crossing  (see  Figure  6).  Typically,  the  time  for 
a  number  of  zero  crossings,  say  8,  is  compared  to  the  time  for  the  first  5  zero 
crossings.  If  these  times  do  not  indicate  the  same  frequency  within  some  allow¬ 
able  error,  e.g.,  2%,  then  the  Doppler  burst  is  rejected  for  further  processing. 

A  block  diagram  of  the  dual-counter  Doppler  signal  analysis  system  is  shown 
in  Figure  7.  The  input  stage  is  a  bank  of  analog  filters.  The  band-pass  and 
rolloff  of  the  analog  filters  were  chosen  to  assure  the  pedestal  would  be 
removed.  The  output  of  each  filter  feeds  a  dual-counter  section  where  each 
counter  counts  a  high-frequency  clock  gated  by  a  selectable  number  of  Doppler 
counts  (3,  5,  8,  or  13  counts).  The  outputs  of  the  dual  counters  are  fed  to 
the  microprocessor,  which  reads  each  of  the  dual-counter  channels  and  deter¬ 
mines  the  highest  frequency  Doppler  channel  that  has  valid  data.  A  valid 
channel  is  one  that  has  had  at  least  8  counts  (if  strapped  for  a  5/8  compari¬ 
son)  and  that  has  a  count  ratio  between  the  two  counters  within  a  software 
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Figure  ft.  Counter-Period  Timing  System  for  Doppler  Signal  Analysis  -  Upper  and 
lower  thresholds  are  used  to  identify  valid  zero  crossing.  A  high- 
frequency  clock  is  counted  for  N  valid  zero  crossings  (gate  time). 
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of  Doppler  Signal  Processing  Electronics  for  Dual  Counter 
am  shown  for  one  velocity  axis.  Additional  analog  circuits 
sed  for  further  axes  with  the  same  CPU  and  operator  interface 


selectable  percentage  of  the  5/8  ratio.  All  of  the  Doppler  channels  are 
sampled  for  valid  data  and  sent  via  IEEE-488  at  the  sample  frequency  selected 
by  the  operator  via  the  host  computer  (COMPAQ). 

System  design 

The  first  step  was  to  produce  an  overall  system  design  for  an  HRC  including 
a  three-component  DLDV  and  temperature  and  conductivity  sensors.  A  detailed 
block  diagram  of  the  system  was  completed  that  identified  off-the-shelf 
components  that  could  be  used.  The  data  processing  electronics  was  based  on 
the  industry  standard  STD  bus  (8  bit).  This  allowed  the  use  of  off-the-shelf 
microprocessor  boards,  remote  computer  interfaces  and  interfaces  to  the 
temperature  and  conductivity  sensor  digital  inputs.  To  retain  compatibility 
with  flow's  microprocessor  development  system  and  the  STD  bus,  an  8088-based 
microprocessor  was  used. 


Prototype  hardware  design 

This  is  part  of  the  circuitry  that  separates  the  pedestal  from  the  Doppler 
signal.  Computer  simulation  was  used  to  select  the  required  bandwidth  and 
rolloff  for  the  analog  front  end.  The  computer  was  used  to  generate  a  worst 
case  Doppler  burst,  and  the  simulation  was  used  to  select  the  bandwidth  and 
rolloff  necessary  to  remove  the  pedestal  and  still  keep  circuit  complexity  to 
a  minimum.  A  rolloff  of  24  dB/octave  on  both  ends  of  the  band-pass  is  suffi¬ 
cient  to  remove  the  pedestal.  The  frequency  ranges  (band-pass)  selected  are 
shown  in  Table  1.  There  are  six  counter  boards  per  cluster  (two  per  axis) 


Table  1.  Frequency  Ranges  for  Counter  Channels 


Channel 

1 

2 

3 

4 

5 

6* 


Frequency  Range 
250.0  Hz  to  1.0  kHz 
1.0  kHz  to  4.0  kHz 
4.0  kHz  to  16.0  kHz 
16.0  kHz  to  64.0  kHz 
64.0  kHz  to  250.0  kHz 
250.0  kHz  to  1.0  MHz 


♦Not  used  for  this  experiment. 
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with  three  dual-counter  channels  per  board.  This  gives  six  frequency  ranges 
per  axis.  A  diagram  of  the  counter  circuit  is  given  in  Appendix  B. 

Each  counter  channel  was  strapped  to  generate  a  5/8  count  ratio.  The  count 
clock  is  set  for  approximately  20  times  the  maximum  frequency  for  a  given  range. 
This  gives  plenty  of  frequency  resolution  over  an  8  Doppler  count  window.  Each 
of  the  counter  boards  fits  into  the  STD  bus  card  cage,  along  with  the  micro¬ 
processor  boards  and  the  temperature  and  conductivity  interface. 

Each  set  of  two  counter  boards  is  driven  by  a  board  mounted  outside  the  STD 
bus  rack  containing  the  analog  signal  processing  electronics  (there  are  three 
analog  boards  per  cluster).  This  board  consists  of  a  differential  input  stage 
driven  by  the  velocimeter  head  for  that  particular  axis.  This  is  followed  by 
the  six  sets  of  filter  stages  (one  for  each  of  the  six  velocity  ranges)  and  six 
comparator  and  gating  stages,  which  convert  the  analog  Doppler  bursts  into 
digital  bursts.  These  digital  bursts  are  then  sent  differentially  to  the  two 
counter  cards  in  the  STD  bus.  There  are  two  adjustments  on  each  of  the  six 
channels:  one  sets  the  amplification  of  the  analog  signal  processing 
electronics,  and  one  sets  the  comparator  threshold  level.  A  9-volt  regulator 
is  also  on  each  analog  board  for  driving  the  laser  diode  (see  Appendix  B). 

Prototype  software  design 

Simultaneously  with  the  hardware  design,  a  software  package  was  developed. 
The  majority  of  the  software  was  written  in  C,  with  a  little  in  assembly 
language.  This  speeds  the  development  time  and  makes  the  software  easier  to 
follow  for  documentation  purposes  and  easier  to  upgrade  for  future 
modifications.  The  software  performs  the  tasks  of  reading,  validating,  and 
interpreting  the  raw  Doppler  signal  data  from  the  dual  counters.  The  software 
is  capable  of  ave -aging  over  several  readings,  keeping  track  of  particle  count 
rates,  and  transmitting  the  data  to  an  auxiliary  processor/microcomputer  for 
data  storage  and  analysis. 


The  unit  is  a  single  circuit  board  designed  to  be  STD  bus  compatible.  The 
prototype  board  performs  three  basic  functions: 

(1)  To  convert  temperature  to  a  voltage  acceptable  to  the  analog  interface 
board  . 
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(2)  To  convert  the  conductivity  of  ocean  water  to  a  voltage  acceptable  to 
the  analog  interface  board. 

(3)  To  provide  prograiranable  control  of  both  temperature  and  conductivity  to 
enable  computer  control  of  gain  and  offset  for  both  temperature  and 
conductivity.  This  capability  provides  very  high  resolution  over  a 
wide  dynamic  range  . 

A  detailed  design  of  the  circuit  board  is  given  in  Appendix  C.  In  essence, 
the  temperature  circuit  is  a  simple,  balanced  resistance  divider  where  the 
temperature-sensitive  element  is  connected  to  a  precision  load  resistor  (Rl). 
Voltage  to  this  divider  is  provided  by  a  precision  voltage  reference  source. 

The  divider  resistor  values  have  been  chosen  to  provide  maximum  sensitivity  at 
the  anticipated  temperature  range  of  -4  to  +15®C.  The  balanced  temperature 
sensor  input  is  preamplified  through  a  very  low-noise,  low-drift  precision 
difference  voltage  follower.  The  signal  is  then  low-pass  filtered  at  either 
3.8,  6.3,  or  12.5  Hz.  The  signal  is  then  a  voltage  representing  temperature 
and  is  controllable  by  the  computer  through  a  selectable  gain  amplifier  and 
off-settable  by  a  16-bit  digital-to-voltage  converter.  The  possible  gains  are 
xl ,  xlO,  xlOO,  and  xlOOO.  Prior  to  the  gain  control,  the  signal  is 
off-settable  from  0  VDC  to  +10.0000  VDC  in  increments  of  0.000153  VDC. 

Conductivity  is  determined  by  the  root-mean-square  (rms)  value  of  an 
alternating  current  (Isig)  flowing  through  the  sample  solution  being  tested. 

The  conductivity  is  equal  to  the  probe  current  (Isig)  divided  by  the  probe 
voltage  (Vsig).  The  circuit  measures  both  Isig  and  Vsig  and  divides  Isig  by 
Vsig  to  obtain  a  voltage  proportional  to  the  conductivity  independent  of 
fluctuations  in  Vsig  (within  limits).  This  voltage  is  then  processed  through 
exactly  the  same  circuitry  as  the  temperature  so  that  the  control  circuitry  has 

t 

a  matched  frequency  response  and  phase  shift. 

Note  that  the  rms  DC  converters  in  the  conductivity  circuit  may  produce  a 
small  correlation  error  of  the  temperature  versus  the  conductivity.  They  also 
limit  the  frequency  response,  but  compared  to  the  low-pass  filters  should  be 
negligible . 


20 


TR-352/11-85 


2,3.3  Miscellaneous  Hardware 


The  underwater  HRCs  are  connected  to  the  processor  through  an  underwater 
cable.  At  the  lower  end  of  the  cable,  a  set  of  pigtails  is  potted  to  a 
junction  box.  At  the  end  of  the  pigtails,  different  types  of  underwater 
connectors  are  plugged  into  individual  sensors  on  the  HRC.  A  total  of  three 
underwater  cables  are  being  fabricated.  These  cables  are  15,  20  and  65  m  long, 
respectively.  The  different  cable  lengths  will  provide  flexibility  for 
measurements  at  different  depths  under  an  ice  floe. 


3.  PERFORMANCE  CHARACTERISTICS  OF  THE  HRC  SENSORS 

Several  series  of  laboratory  tests  were  conducted  during  Phase  I  and 
Phase  II  to  determine  the  performance  characteristics  of  the  sensors  used  on 
the  HRCs.  The  tests  conducted  in  Phase  I  are  reported  in  detail  in  [I].  In 
this  section,  we  describe  the  Phase  II  tests  and  results.  From  these  results, 
we  define  the  performance  characteristics  of  the  sensors. 


3.1  Diode  Laser  Doppler  Velocimeter 

3.1.1  Prototype  DLDV 

For  a  dual-beam  system  using  the  fringe  model  [I],  the  calibration  of  the 
LDV  depends  on  the  wavelength  of  the  laser,  X,  and  the  half  angle  of  the  beam 
crossing,  <|),  or 


2  u  sin  4» 

X 


(1) 


where  f^  is  the  Doppler  frequency  and  u^  is  the  velocity  component 
perpendicular  to  the  bisector  of  the  two  illuminating  beams.  Replacing 
X/(2  sin  <j>)  =  dj,  we  have 


[Hz/ (cm/ s) ] 


(2) 


where  d^  is  the  fringe  spacing. 

These  two  equations  indicate  that  the  calibration  of  the  DLDV  is  linear 
and  its  sensitivity  in  hertz  per  centiroeters/second  is  inversely  proportional 
to  the  fringe  spacing.  For  a  DLDV  with  a  known  geometrical  arrangement  of  the 
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crossing  beams,  the  calibration  will  stay  the  same  until  the  laser  diode  is 
replaced.  For  a  given  diode,  the  peak  wavelength  (about  790  nm  for  a  Hitachi 
HL7801G  diode)  increases  slightly  with  the  output  power  and,  therefore,  with 
the  device  (but  not  the  ambient)  temperature.*  The  typical  average  temperature 
coefficient  of  the  peak  wavelength  is  0.2  to  0.3  nm/°C.  In  the  extreme  situa¬ 
tion  of  a  10®C  drop  in  the  device  temperature  in  response  to  a  20'’C  ambient 
temperature  drop,  the  peak  wavelength  decreases  less  than  0.4%.  This  corre¬ 
sponds  to  an  increase  of  0.4%  in  the  sensitivity,  which  may  be  corrected  for 
if  the  device  temperature  is  known.  In  practice,  such  a  small  correction  is 
generally  well  within  the  experimental  error  (typically  +  5  to  10%)  for  field 
experiments . 

A  series  of  tow-tank  experiments  was  conducted  at  FLOW  prior  to  MIZEX  84 

to  confirm  Equation  (2).  One  of  the  DLDV  units  tested  had  an  HL7801G  diode 

(No.  0005)  with  a  peak  wavelength  of  793  nm  at  3  mW.  The  half  angle  measured 

in  water,  (j)  ,  was  3.5°.  The  fringe  spacing  was  therefore 
w 

d  =  X  /2  sin  (j)  =  X  /2  n  sin  d)  =  4.87  (3) 

fa  ^a  a  w  ^w 

where  n^  =  1.333  is  the  index  of  refraction  of  water.  The  corresponding 
sensitivity,  according  to  Equation  (2),  was  estimated  to  be  2052  Hz/(cm/s). 
Figure  8  shows  the  results  of  the  calibration.  The  ordinate  and  abscissa  are, 
respectively,  the  velocity  of  the  tow  carriage  on  which  the  DLDV  was  mounted, 

U_,  and  that  measured  with  the  DLDV,  U  ,  As  can  be  seen  in  the  figure, 
the  data  points  are  within  ^  3%  of  the  45°  line.  Similar  results  were  obtained 
with  another  DLDV  unit  in  which  an  LDL  SCW21  diode  (831  nm  at  7  mW)  was 
installed.  Both  units  were  used  during  MIZEX  84. 

Additional  tow-tank  tests  were  conducted  to  test  the  three  3-axis  clusters 
prior  to  AIWEX  85.  For  all  nine  units,  the  nominal  wavelength  of  the  diodes 
was  790  ^  5  nm,  and  the  nominal  half  angle  between  the  beams  was  6.71°  in 
air.  The  fringe  spacing  was  therefore  3.38  v»n  and  the  sensitivity  was 
2960  Hz/ (cm/s) . 

During  the  tow  tests,  we  had  to  add  powdered  milk  into  the  tank  in  which 
the  water  was  filtered  through  a  lO-^m  filter.  There  were  simply  not  enough 
particles  left  in  the  filtered  water  for  the  DLDV  to  work  satisfactorily . 

♦Sources  quoted  from  "Optical  Semiconductor  Devices  Data  Book"  issued  by 
Hitachi  America,  Ltd.,  Documentation  No.  S13. 


22 


TR-352/11-85 


(cm/s) 


3.1.2  Configuration  for  Three-Axis  HRC  System 

As  shown  in  Figure  1,  three  single-axis  DLDVs  are  mounted  on  an  octagonal 
Itame  that  is  held  vertically.  The  transmitter  receiver  pairs  are  identified 
as  #1,  #2  and  #3.  U1  measures  the  velocity  component  perpendicular  to  the 
plane  of  the  frame,  while  U2  and  U3  measure  the  velocity  components  in  the 
plane  of  the  frame  inclined  at  45°  to  the  vertical.  The  frame  is  rotated  in 
the  horizontal  plane  so  that  the  current  flows  into  the  frame  at  an  angle  of 
30  to  45®  to  the  right  of  the  frame  axis.  This  ensures  that  all  three  DLDV 
pairs  observe  a  substantial  mean  flow  component.  The  orientation  of  all  three 
DLDV  pairs  is  shown  in  Figure  1. 

3 . 2  Temperature  and  Conductivity  Probes 

The  FASTIP  thermistor  by  Thermometries  (Series  FP07)  has  a  very  thin  tip 
(0.18  mm  overall  diameter)  in  which  the  thermistor  bead  is  housed.  As  a 
result,  it  has  a  frequency  response  of  21  Hz,  which  is  about  a  factor  of  2 
higher  than  those  sealed  in  a  glass  tube  without  a  thin  tip  (Lueck  et  al . , 
1978;  Gregg  and  Meagher,  1980).  The  thermistor  circuit  was  designed  to  have 
the  maximum  sensitivity  at  temperatures  around  0°C  (see  Appendix  C) .  At  a 
gain  of  1,  the  sensitivity  was  typically  5  to  6  mV/®C.  Due  to  the  relatively 
low  sensitivity,  any  small  electronic  drift  (not  the  sensitivity)  would  result 
in  a  noticeable  shift  in  the  measured  absolute  temperature.  Therefore,  this 
particular  design  was  not  suitable  for  absolute  temperature  measurement. 
However,  the  absolute  temperature  may  be  inferred  from  nearby  conventional 
temperature  gauges  deployed  simultaneously.  The  rms  noise  level,  which  was 
predominantly  60  Hz,  was  measured  to  be  about  0.7  mV  at  a  gain  of  100  and  with 
the  low-pass  filter  set  at  100  Hz.  The  noise  level  decreased  significantly 
when  the  low-pass  filter  was  subsequently  set  at  or  below  12.5  Hz.  Therefore, 
the  thermistor  gauge  is  capable  of  resolving  minute  temperature  changes  of 
about  0.001®C.  The  sensitivity  of  the  thermistor  varies  with  the  temperature 
but  may  be  considered  constant  within  a  narrow  temperature  range  of  2  to  3°C  . 
The  actual  resolution  of  the  gauge  is  about  0.002®C  at  a  gain  of  100  and  a 
12-bit  A/D  conversion  (10  V  full  scale). 

As  described  in  Section  2.2,  the  dual-electrode  conductivity  probes  have 
been  subjected  to  field  tests  and  found  to  be  satisfactory  for  measurements  of 
conductivity  differentials  (Meagher  et  al.,  1982).  To  ensure  sufficient 
spatial  resolution  (1  to  2  cm),  we  reduced  the  diameter  of  the  platinum  wires 
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2.5  times  from  1.3  mm  to  0.5  mm  at  the  sacrifice  of  some  sensitivity  and, 
probably,  stability.  The  frequency  response  of  the  dual-electrode  probe  is 
expected  to  reach  hundreds  of  hertz  (Dugan  and  Okawa,  1982).  To  avoid  salinity 
"spiking"  (Dantzler,  1974;  Roden  and  Irish,  1975),  however,  the  circuit  has  the 
same  low-pass  filter  as  that  of  the  thermistor  probe,  i.e.,  12.5  Hz  or  lower. 

The  circuitry  of  the  dual-electrode  probe  was  designed  to  have  voltage 
output  linearly  proportional  to  conductivity.  Calibration  of  the  conductivity 
gauge  was  made  by  using  four  jars  of  salt  solution  with  different  brine 
concentrations  (p  =  1.0,  1.01,  1.02,  and  1.025  gr/cm  ).  The  jars  were  kept 
in  a  constant  temperature  bath  at  20“C.  The  conductivity  was  measured  with  a 
Beckman  conductivity  gauge  (Model  RC-16C)  using  a  standard  cell  with  a  cell 
constant  of  1.0  cm  Figure  9  shows  three  repeated  sets  of  calibration 
data.  The  nominal  sensitivity  is  about  15  mmho/cm/V  at  a  gain  of  1.  This 
sensitivity  deteriorates  slowly  with  time  in  submersion  due  to  fouling  of  the 
platinum  black  coating.  Typically,  the  sensitivity  decreases  2  to  5^  after  2 
to  3  days  of  submersion  in  the  calibration  jars.  To  revive  a  desensitized 
probe,  one  may  simply  dip  it  in  a  10%  HCl  solution  for  2  to  3  minutes  (a 
standard  practice  recommended  by  Beckman  for  their  conductivity  cells).  One 
may  also  replatinized  the  probe,  but  that  would  change  the  sensitivity 
noticeably  depending  on  the  condition  of  the  replatinized  tip  surfaces,  which 
cannot  be  controlled  accurately. 

At  a  gain  of  100  and  with  the  low-pass  filter  set  at  100  Hz,  the  rms  noise 
level  was  measured  to  be  about  26  mV.  This  noise  level  again  decreased  notice¬ 
ably  after  the  low-pass  filter  was  set  at  or  below  12,5  Hz.  At  a  gain  of  100, 
the  dual-conductivity  probe  is  capable  of  resolving  minute  conductivity  changes 
of  1  V»nho/cra.  The  actual  resolution  is  about  2  yrnho/cm  at  a  gain  of  100 
and  a  12-bit  A/D  conversion  (10  V  full  scale). 


4.  SOFTWARE  FOR  DATA  ACQUISITION  AMD  ANALYSIS 
4.1  Data  Acquisition  Software 

The  V-T-C  HRC  is  an  intelligent  instrument  that  includes  a  GPIB  interface 
for  transfer  of  data  to  a  host  computer.  We  designed  the  software  to  maximize 
use  of  the  interface  so  that  all  system  parameters  were  under  direct  program 
control.  Temperature  and  conductivity  probe  calibration  coefficients  are 
downloaded  to  the  cluster  prior  to  the  experiment.  The  user,  through  the  host 
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Figure  9.  Calibration  Data  for  Conductivity  Probes 


computer  cluster  communications  program,  may  also  specify  a  sampling  rate  from 

3  to  25  Hz.  When  a  trigger  command  is  sent  to  the  HRC,  sampling  of  the  three 
DLDV  velocity  channels  and  the  high-resolution  temperature  and  conductivity 
channels  begins.  After  each  sampling  time,  the  values  of  the  three  velocity 
components,  temperature  and  conductivity  are  transferred  over  the  GPIB  to  the 
host  computer  for  storage  on  a  hard  disk. 

Several  independent  clusters  may  be  connected  in  parallel  on  the  GPIB  bus, 
so  that  a  single  host  computer  may  simultaneously  accept  data  from  all 
clusters.  The  computer  software  was  designed  to  be  very  flexible  in  order  to 
allow  users  to  configure  the  system  to  their  requirements. 

The  system  is  controlled  by  an  R088-based  microprocessor  on  an  STD  bus 
card.  An  8087  numeric  data  processor  also  resides  on  the  card  to  carry  out 
fast  floating-point  calibration  operations.  The  operating  instructions  for 
the  data  acquisition  software  are  given  in  Appendix  A. 

4  .2  Data  Analysis  Software 

The  raw  data  files  were  too  large  to  store  on  floppy  disk,  so  a  program 
called  CULL  was  written  to  read  them  and  do  some  primary  processing.  The  raw 
data  from  different  processors  were  stored  in  an  interleaved  fashion  on  a 
single  file.  The  CULL  program  first  separates  the  data  according  to  cluster 
and  then  determines  the  time  at  which  each  data  point  was  taken.  Null  data 
are  eliminated,  and  Doppler  frequencies  are  transformed  into  velocities. 

A  number  of  problems  were  encountered  with  spurious  data  points  in  the 
files  resulting  from  inadequate  filtering  and  discrimination  by  the  processor. 
The  CULL  routine  thus  has  a  provision  for  a  minimum  acceptable  velocity  from 
each  channel.  In  addition,  a  persistence  check  is  made.  Sequential  data 
points  are  assumed  to  be  invalid  if  they  are  outside  a  fixed  range  about  the 
last  valid  data  point.  This  persistence  was  generally  set  at  2  to  3  times  the 
variance  of  the  velocity.  As  time  went  on  (during  AIWEX  85)  and  the 
processors  were  better  adjusted  and  the  proper  filterbands  disabled,  the 
problem  with  spurious  data  points  abated. 

The  files  created  by  CULL  were  on  the  order  of  300  kilobytes  long  so  that 
they  could  be  stored  on  floppy  disk.  These  files  were  also  short  enough  that 
it  was  possible  to  read  through  them  using  a  text  editor  to  check  for  spurious 
data,  data  rate  and  overall  system  behavior.  Some  of  the  files  were  too  long 
to  be  edited  and  were  lost  when  this  was  attempted  . 
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The  CULL  routine  counts  the  number  of  valid  data  points  on  each  channel  and 
writes  this  information  at  the  top  of  the  output  file.  In  addition,  it  was 
necessary  to  reinput  the  temperature/conductivity  probe  parameters  because  the 
binary  files  were  written  with  tabs,  which  are  not  readable  in  Fortran.  CULL 
requires  approximately  one  hour  of  processing  time  on  an  average-sized  file. 

Two  programs  were  written  to  provide  initial  analysis  of  the  data.  These 
provide  Reynolds  stress,  turbulent  fluxes  of  heat  and  salinity,  variances  of 
velocity,  temperature  and  conductivity  and  the  one-dimensional  wavenumber 
spectra  of  velocity  fluctuations. 

The  first  program  is  called  TURB.  This  program  separates  the  velocity 
into  a  mean  and  a  fluctuating  part, 

I 

u.  =  U .  +  u  . ,  (4 ) 

111 


and  then  calculates  the  variance-covariance  tensor  of  the  three-component 
fluctuating  velocity  field. 
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The  overbar  indicates  averaging  over  the  observed  values  of  velocity.  The 
diagonal  elements  of  the  tensor  are  the  variances  of  the  three  observed  velo¬ 
city  components,  while  the  off-diagonal  terms  are  covariances.  The  covariance 
terms  describe  the  shear  stress  in  the  flow.  These  must  be  rotated  into  a 
reference  frame  with  one  axis  vertical,  one  in  the  horizontal  plane  aligned 
with  the  mean  flow  and  the  last  orthogonal  to  these  two.  This  is  accomplished 
with  two  rotation  matrices.  The  first  provides  a  reference  frame  rotation  of 
45®  in  a  counterclockwise  direction  about  the  horizontal  axis  of  the  cluster 
so  that  the  U3  coordinate  axis  is  vertical  and  pointing  down.  This  rotation 
matrix  is 
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The  next  rotation  is  about  the  vertical  axis  to  provide  a  reference  frame 
aligned  with  the  mean  flow  direction  in  the  horizontal  plane: 
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The  complete  rotation  of  the  variance-covariance  tensor  is  accomplished  by 


a’  = 


^  V2 


(9] 


This  tensor  gives  the  Reynolds  stress  for  vertical  momentum  flux  or,  equiva¬ 
lently,  the  friction  velocity: 
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(10) 


The  same  program  obtains  the  turbulent  flux  vectors  of  heat  and  salinity. 


and 


=  h'u‘ 
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The  fluctuating  component  of  sensible  heat  is  obtained  from  the  temperature 
fluctuations 


H'  =  c  T'  (13 

P 

where  c^  is  the  heat  capacity  of  water.  The  fluctuating  part  of  salinity 
was  obtained  by  numerically  differentiating  the  relationship  given  in  "1978 
Practical  Salinity  Scale  Equations,"  IEEE  Journal  of  Ocean  Engineering, 

Vol .  OE-5,  No.  1,  January,  1980.  This  algorithm  was  used  to  provide  values  of 
dS/dC  at  constant  temperature,  pressure  and  absolute  conductivity,  resulting  in 
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and  are  vectors  along  the  axes  of  the  DLDV  cluster  reference  frame. 
Rotation  of  these  vectors  to  obtain  vertical  fluxes  is  accomplished  with  the 
same  matrices  used  to  rotate  the  variance-covariance  tensor, 


x'  = 


(15) 


The  variance  of  the  three  components  of  motion  as  well  as  temperature  and 
conductivity  is  also  obtained.  This  program  requires  over  300  kilobytes  of 
random  access  memory  to  accommodate  the  large  data  arrays.  Processing  time  is 
about  10  minutes  per  cluster. 

A  program  was  also  written  to  determine  the  wavenumber  spectrum  of  each 
velocity  component  as  well  as  temperature  or  conductivity.  The  time  series  of 
a  given  record  is  transformed  into  the  spatial  dimension  by  multiplying  by  the 
mean  flow  velocity.  It  is  thus  assumed  that  the  turbulence  is  horizontally 
isotropic  and  that  the  statistics  are  stationary  over  the  period  of  the 
record.  This  is  known  as  the  Taylor  frozen  turbulence  hypothesis.  Laser 
Doppler  velocimeter  data  are  characterized  by  time  series  of  velocities 
obtained  at  random  time  intervals.  This  makes  the  spectral  analysis  using  the 
standard  FFT  approach  impossible  because  the  FFT  algorithm  requires  that  the 
time  series  to  be  transformed  must  be  obtained  at  regular  intervals.  An 
alternative  technique  was,  therefore,  used  for  the  DLDV  data.  The  spectrum  is 
calculated  using  the  Wiener-Khinchine  theorem,  which  states  that  the  spectral 
density  is  given  hy  the  Fourier  transform  of  the  autocovariance  function.  The 
autocovariance  function  is  obtained  from 


with 


R. .  =  u(x. )  u(x. +r) 
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where  r  is  the  wavenumber  offset  analogous  to  a  time  lag  in  frequency  space. 
This  function  can  be  estimated  for  time  series  that  are  spaced  at  random  time 
intervals  by  dividing  the  time  lag  domain  into  discrete  slots.  The  estimator 
is 
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where  H(r)  is  the  number  of  lagged  products  obtained  when  the  time  lag  lies 
within  the  time  slot,  r.  The  time  slot  for  a  given  product  is  defined  by 

<0.5  .  (18) 

The  averaging  is  done  numerically  using  the  slotted  correlogram  technique 
described  by  Bell  (1983).  Bell  compared  this  technique  with  other  approaches 
to  spectral  analysis  of  laser  velocimetry  and  found  it  to  be  the  best  compro¬ 
mise  between  accuracy  and  efficiency.  The  spectrum  is  then  calculated  using 
the  FFT.  Spectral  densities  are  normalized  by  the  variance  of  the  motion  to 
allow  direct  intercomparison  and  averaging  of  spectra.  Errors  in  the  spectrum 
are  related  to  the  number  of  cross  products,  H,  utilized  in  the  average  value 
at  each  offset  value.  The  signal-to-noise  ratio  of  the  spectral  density  at 
wavenumber  k  is 

SNR  =  yiilT/rk 

The  program  used  to  do  this  is  called  SPFAST.  It  requires  a  RAM  of  over 
300  kilobytes  and  takes  about  an  hour  to  provide  one  spectrum. 

3 .  FIELD  MEASOREMENTS  AND  RESULTS 

To  demonstrate  the  field  deployability  of  the  DLDVs,  we  built  two  experi¬ 
mental  units  and  conducted  a  series  of  tests  during  MIZEX  84.  Documentation 
of  the  MIZEX  84  tests  is  given  in  Section  5.1.  From  the  results  of  MIZEX  84, 
which  were  very  encouraging,  we  optimized  the  design  of  the  DLDV  modules  and 
fabricated  nine  sets  of  prototype  DLDV  modules  for  use  in  three  3-axis  high- 
resolution  clusters.  In  addition  to  the  three  DLDVs,  each  cluster  included  a 
high-frequency  thermistor  and  a  conductivity  probe.  Each  HRC  was  controlled 
by  a  processor  that  provided  power  to  drive  the  sensors,  preprocess  the  raw 
data  and  interface  the  sensor  output  to  the  host  computer.  The  three  HRCs 
were  then  deployed  during  AIWEX  85.  Documentation  of  AIWEX  85  tests  is  given 
in  Section  5.2. 
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5.1  MIZEX  84 


Two  field-deployable  DLDVs  were  assembled  and  field  tested  during  MIZEX  84. 
The  objectives  of  our  participation  were  twofold:  (1)  to  determine  the  feasi¬ 
bility  of  the  r»LDV  for  use  in  the  Arctic  environment,  which  is  the  most  severe 
test  to  the  field  deployability  of  this  optical  sensor,  and  (2)  to  gather 
design  information  for  the  prototype  HRCs  using  the  DLDVs  as  the  primary 
velocity  sensors,  provided  the  feasibility  was  successfully  demonstrated. 

5.1.1  MIZEX  84  Site  and  Setup 

The  MIZEX  field  site  was  located  well  within  the  marginal  ice  zone  at  about 
80®  north  latitude  and  0®  west  longitude  in  the  Frara  Straits  of  the  Greenland 
Sea.  Living  quarters  and  workshops  were  provided  on-board  the  M/V  Polar  Queen. 
The  DLDV  tests  were  carried  out  on  a  2-  to  3-m-thick  floe  at  a  distance  of 
about  250  m  from  the  ship.  The  moderate-resolution  turbulence  measurement 
experiment  carried  out  by  Miles  McPhee  was  set  up  at  the  same  location. 

A  1.2-m  square  hole  was  prepared  in  the  ice  for  the  DLDV  frame.  The  ice 
at  this  point  was  2.4  m  thick.  The  electronic  equipment  was  set  up  in  a 
heated  tent  nearby.  Two  DLDV  systems  were  deployed  during  the  test.  Power 
was  provided  by  a  600-W  gasoline-powered  generator.  The  DLDV  systems  required 
two  DC  power  supplies.  Outputs  from  the  photodetectors  were  filtered  through 
two  Krohn-Hite  variable  band-pass  filters.  The  filtered  Doppler  signals  were 
processed  into  analog  signals  proportional  to  the  current  by  two  systems,  a 
Macrodyne  Model  2096  Laser  Doppler  Signal  Processor  and  a  Thermal  Systems  Inc. 
(TSI)  Model  1090  Laser  Doppler  Frequency  Tracker,  that  used  different  pro¬ 
cessing  schemes.  The  analog  signal  was  sampled  at  a  rate  of  48  Hz  and 
converted  into  a  digital  signal  by  a  modified  Sea-Bird  Electronics  Deck  Unit. 
The  digital  data  were  recorded  onto  3.5-inch  floppy  disks  using  a  Hewlett- 
Packard  9816  microcomputer.  Doppler  signals  were  observed  directly  using  a 
dual-channel,  Tektronics  storage  scope. 

The  DLDV  laser  units  were  aligned  at  0®C  in  a  shack  built  on  the  ice. 

This  is  necessary  to  avoid  thermal  strains  in  the  optics,  which  result  in  the 
beams  no  longer  crossing  and  the  Doppler  signal  being  lost  when  the  system  is 
lowered  into  cold  water.  Only  one  of  the  systems  was  found  to  drift  out  of 
alignment  when  cooled  from  25  to  0®C.  This  system  also  went  out  of  alignment 
when  brought  to  room  temperature  and  returned  to  0®C. 
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The  laser  units  and  photodetectors  were  then  mounted  on  a  rigid  tubular 
frame  and  suspended  just  beneath  the  water  surface  for  alignment  with  each 
other.  Underwater  alignment  was  necessary  because  the  front  surface  of  the 
focusing  lens  on  the  laser  unit  is  nonplanar.  Since  the  curvature  of  this 
lens  was  unknown,  the  relationship  between  the  beam  crossing  angle  in  air  and 
underwater  was  unknown.  The  beam  crossing  angle  is  used  to  determine  the 
Doppler  fringe  spacing  and,  hence,  the  relationship  between  Doppler  frequency 
and  measured  current  speed.  The  beam  crossing  angles  were  observed  to  be  the 
same  as  measured  in  the  laboratory.  The  use  of  a  lens  with  a  planar  front 
surface  is  recommended  in  the  future  because  this  allows  the  fringe  spacing  to 
be  calculated  from  measurements  in  air  and  simplifies  the  locating  of  the 
crossing  point  (Doppler  fringe  volume)  underwater. 

A  fiber  optics  rod  3  mm  in  diameter  is  used  to  transmit  light  from  the 
measurement  point  to  a  photodetector  diode.  This  rod  must  be  located  so  that 
it  transmits  the  light  that  is  forward  scattered  by  particles  passing  through 
the  sampling  volume.  The  alignment  must  be  carried  out  in  water  at  0®C.  This 
was  done  by  suspending  the  DLDV  frame  just  underneath  the  water  in  the 
hydro-hole.  The  experiment  was  carried  out  during  the  height  of  the  Arctic 
summer,  so  it  was  not  possible  to  observe  the  laser  beams  outdoors.  A  shack 
was  therefore  constructed  over  the  hydro-hole  to  allow  positioning  of  the 
photodetector  fiber  optic  rod.  The  end  of  the  rod  was  located  6  cm  from  the 
crossing  point  and  centered  between  the  two  crossed  laser  beams.  The  laser 
units  and  photodetectors  were  mounted  on  adjustable  mounts  that  allowed  for  a 
straight-forward  alignment.  Once  the  photodetector  appeared  to  be  properly 
aligned,  a  check  was  made  to  see  if  a  Doppler  signal  could  be  obtained.  This 
was  done  by  stirring  the  water  in  the  hole  while  observing  the  photodetector 
output  on  a  storage  oscilloscope.  Better  alignment  could  be  made  in  a  tank  of 
ice  water,  so  that  the  Doppler  signal  could  be  observed  while  the  alignment 
was  carried  out.  In  the  field,  this  would  require  a  sizeable  instrument  shack 
capable  of  being  darkened  to  allow  observation  of  the  laser  light  and  having 
enough  space  for  all  the  electronics.  During  the  Arctic  winter,  an  enclosed, 
heated  space  will  be  required  for  field  alignment. 

Once  good  Doppler  signals  were  obtained  from  both  DLDV  systems,  the  systems 
were  lowered  beneath  the  ice.  The  ice  thickness  in  the  hole  was  2.4  m,  and  the 
array  was  lowered  2  m  below  the  bottom  surface  of  the  ice.  This  is  a  standard 
depth  for  under-ice  measurements  of  current,  and  is  generally  well  within  the 
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boundary  layer.  One  of  the  systems  was  oriented  to  observe  the  horizontal 
flow,  while  the  other  was  oriented  to  observe  the  component  of  flow  at  45®  from 
the  vertical.  A  tell-tale  was  used  to  orient  the  frame  into  the  mean  flow  so 
that  any  turbulent  fluctuations  would  be  superimposed  onto  the  mean.  This  is 
done  to  resolve  the  directional  ambiguity  in  the  observed  current.  The  cali¬ 
bration  constants  and  orientations  of  the  two  DLDVs  are  given  in  Table  2. 


Table  2.  Calibration  Constants  and  Orientations  of  DLDVs 

Laser 
Hitachi 
LDL 

The  output  power  from  a  laser  diode  increases  as  the  operating  temperature 
is  reduced.  If  the  output  power  exceeds  a  critical  level,  the  laser  will  be 
damaged.  The  Hitachi  laser  diode  is  equipped  with  a  photodiode  monitor  that 
was  used  in  a  feedback  loop  to  maintain  constant  laser  power  output  regardless 
of  temperature.  The  laser  unit  containing  the  Hitachi  laser  was  aligned  at 
room  temperature  and  was  found  to  maintain  its  alignment  after  transport  to  the 
field  site  and  at  freezing  temperatures.  No  field  adjustments  were  necessary 
to  obtain  a  good  Doppler  signal  from  this  unit.  The  LDL  laser  diode  has  no 
monitoring  provision.  The  forward  current  on  this  laser  diode  was  thus  set  at 
a  constant  value  to  give  a  power  output  of  5  mW  at  0®C.  The  laser  unit  con¬ 
taining  the  LDL  laser  diode  was  found  to  be  out  of  alignment  after  transporta¬ 
tion  to  the  field.  This  unit  was  found  to  be  more  sensitive  to  thermal  strains 
so  that  it  was  not  possible  to  align  the  unit  at  room  temperature.  Alignment 
was  therefore  carried  out  in  the  field  shack  at  temperatures  near  freezing. 

Once  this  unit  was  aligned,  a  good  Doppler  signal  was  obtained.  The  strength 
of  the  signal  obtained  from  the  LDL  system  was  never  as  good  as  that  obtained 
from  the  Hitachi  unit,  and  over  a  period  of  a  week  the  Doppler  signal  from  the 
LDL  deteriorated  significantly.  This  deterioration  was  apparently  caused  by  a 
gradual  creep  in  the  glue  joint  of  a  broken  glass  plate  holding  the  beam 
splitter  cubes. 
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A  total  of  2A0  tninutps  of  h i^h-resolut ion  turbulent  flow  data  was  recorded 
onto  floppy  disks,  of  this  total,  42  minutes  were  taken  at  times  when  the 
data  rate  was  hie'  enough  to  obtain  both  the  horizontal  and  45®  components  of 
the  flow  simultane  i^lv.  The  HLOV  units  with  the  LDL  and  Hitachi  diode  lasers 
were  used  for  the  horizontal  and  45®  configurations,  respectively.  Up  to  two 
channels  of  Dl.r",'  data  were  recorded  simultaneously.  The  Macrodyne  counter  was 
always  used  with  channel  0,  whereas  the  TSI  tracker  was  always  used  with 
channel  1.  In  all  of  these  records,  a  simultaneous  recording  was  made  of  the 
output  from  Miles  McPhee's  KRC .  This  array  was  mounted  50  cm  above  the  DLDV 
measurement  volume.  The  magnitude  and  orientation  of  the  mean  horizontal  flow 
was  thus  always  measured.  Data  from  the  two  types  of  current  meters  were  used 
for  comparison  to  assess  their  performance  and  the  importance  of  turbulent 
motions  on  spatial  scales  of  less  than  2  cm. 

For  the  most  part,  data  collection  was  limited  to  one  channel.  Since  the 
mean  horizontal  flow  is  required  to  resolve  the  vertical  component  of  the  45° 
flow,  it  was  usually  not  possible  to  obtain  directly  the  vertical  flow  com¬ 
ponent.  Back-to-back  records  of  horizontal  and  45°  flows  were  taken  to  obtain 
a  reasonable  estimate  of  the  mean  horizontal  flow  for  use  in  resolving  the  45° 
component . 

5.1.2  Results  and  Discussion 

During  MIZEX  84,  a  total  of  18  DLDV  data  records  was  obtained  at  a  depth 
of  2.5  m  below  the  bottom  of  the  ice  floe  and  stored  on  floppy  disks.  Each 
record  contains  14  minutes  of  data  on  two  channels.  Channel  0  is  data  output 
from  the  Macrodyne  counter,  and  channel  1  is  data  from  the  TSI  tracker.  Data 
were  sampled  by  the  computer  at  a  rate  of  48  Hz  for  a  total  of  40,320  data 
points  per  channel.  Much  of  the  data  is  redundant  because  the  Doppler  burst 
rate  was  less  than  48  Hz.  A  look  at  the  raw  data  shows  that  alternate  data 
points  are  slightly  offset  in  time.  This  was  apparently  caused  by  our  attempt 
to  double  the  sampling  rate  available  on  the  Sea-Bird  A/D  converter.  To 
simplify  processing  and  reduce  the  data  files  to  a  manageable  size,  only 
alternate  data  points  were  used  in  the  analysis.  The  effective  sampling  rate 
for  data  analysis  is  thus  24  Hz .  A  number  of  runs  were  made  simultaneously 
with  ducted  propeller  meters  on  an  MRC  placed  50  cm  above  the  DLDV  cluster. 
Table  3  is  a  summary  of  the  DLDV  and  MRC  records. 


Table  3.  Summary  of  DLDV  and  MRC  Records 
(fromMcPhee,  1985b) 


Macrodyne 

(Channel  0) 

TSI 

(Channel  l) 

Run 

Date 

Time 

Sensor 

Orientation 

Sensor 

Orientation 

197  .F 

7/15 

15:26 

LDL 

Hor. 

— 

— 

197. F 

7/15 

15:46 

Hitachi 

45° 

— 

— 

197  .F 

7/15 

16:414 

LDL 

Hor . 

Hitachi 

45° 

198  .A 

7/16 

20:10 

Hitachi 

45° 

— 

— 

198. B 

7/16 

20:34 

LDL 

Hor . 

— 

— 

198. C 

7/16 

22:23 

Hitachi 

45° 

— 

— 

198. D 

7/16 

22:45 

LDL 

Hor . 

— 

— 

199  .A 

7/17 

11:46 

LDL 

Hor . 

199  .B 

7/17 

12:07 

Hitachi 

45° 

— 

— 

199. C 

7/17 

15:37 

Hitachi 

45° 

Hitachi 

45° 

1P9.D 

7/17 

16:00 

LDL 

Hor . 

Hitachi 

45° 

199. E 

7/17 

17:12 

LDL 

Hor . 

Hitachi 

45° 

199  .F 

7/17 

17:33 

Hitachi 

45° 

Hitachi 

45° 

199  .G 

7/17 

20:53 

Hitachi 

45° 

Hitachi 

45° 

199  .H 

7/17 

21:14 

LDL 

Hor . 

Hitachi 

45° 

Doppler  Burst  Rate  and  Data  Rate 

Traces  of  the  Doppler  burst  signals  recorded  on  a  storage  scope,  under  a 
variety  of  conditions,  were  examined  to  understand  the  performance  of  the 
DLDVs  in  the  Arctic  environment.  Doppler  bursts  from  the  photodetectors  were 
band-pass  filtered  to  remove  the  pedestal  and  high-frequency  noise.  For  an 
optimally  aligned  DLDV,  the  typical  signal-to-noise  ratios  (SNRs)  are  better 
than  20:1  after  band-pass  filtering.  At  times,  we  observed  the  presence  of  an 
intermittent  2  MHz  noise  thought  to  be  associated  with  radio  transmissions. 

The  0.5  MHz  low-pass  filter  available  on  the  Macrodyne  counter  reduced  but  did 
not  eliminate  the  noise.  For  a  misaligned  DLDV  with  a  poor  SNR,  this  noise  may 
overwhelm  the  Doppler  signal.  As  a  result,  this  noise  never  caused  problems 
for  the  DLDV  unit  with  the  Hitachi  laser  diode  but  affected  the  signals 
measured  by  the  DLDV  unit  with  the  LDL  laser  diode  when  the  alignment 
deteriorated . 


From  the  outputs  of  the  Macrodyne  counter  and  the  TSI  tracker,  we  estimated 
the  Doppler  burst  rate.  The  most  direct  measure  was  made  by  recording  the 
photodetector  output  signal  using  line  triggering  on  the  storage  oscilloscope. 
The  sweep  rate  was  set  at  2  ms/div,  and  a  total  of  3  seconds  of  signal  was 
examined.  Each  sweep  contained  20  ms  of  randomly  sampled  record.  Doppler 
bursts  are  easily  identified  at  this  sweep  rate.  In  3  seconds,  26  bursts  were 
identified  for  a  burst  rate  of  about  9  per  second.  Immediately  afterwards, 
both  the  Macrodyne  counter  and  the  TSI  frequency  tracker  were  attached  to  the 
photodetector  output,  and  5  seconds  of  output  was  recorded.  The  threshold 
level  of  the  Macrodyne  counter  and  the  gain  on  the  TSI  tracker  were  adjusted 
so  that  both  systems  responded  to  the  same  signal.  This  ensures  that  noise 
signals  are  not  being  interpreted  as  Doppler  bursts.  The  data  rate  taken  from 
this  test  was  about  5  samples  per  second.  The  Macrodyne  system  was  observed 
to  not  trigger  on  Doppler  bursts  that  were  identifiable  as  such  but  were  too 
low  in  magnitude  or  had  too  few  fringes.  The  data  rate  is  thus  substantially 
lower  than  the  Doppler  burst  rate.  The  mean  current  during  this  test  was 
about  7  cm/s,  which  gives  a  spatial  resolution  of  1.4  cm. 

A  50-second  record  of  Doppler  period  was  analyzed.  The  data  rate  varies 
from  3  to  4  per  second.  The  actual  data  rate  is  somewhat  higher  since  the 
oscilloscope  resolution  is  not  fine  enough  to  detect  small  changes  in  the 
output.  The  DLDV  is  oriented  to  measure  the  component  of  flow  45®  from  the 
vertical  in  the  direction  of  the  mean  flow.  The  current  during  this  time 
varied  from  4  to  9  cm/s,  with  a  mean  of  about  6  cm/s;  the  spatial  resolution 
is  thus  about  1.7  cm. 

It  was  not  possible  to  get  consistent  tracking  of  the  Doppler  signal  with 
the  TSI  processor  at  these  low  data  rates.  In  two  of  the  records,  the  same 
DLDV  Doppler  signal  was  fed  into  both  signal  processing  systems.  The  output 
was  converted  into  physical  units  to  give  a  direct  comparison  of  the  two 
processing  techniques.  The  TSI  output  is  offset  by  about  1.5  to  2.5  cm/s  from 
the  Macrodyne  signal  and  exhibits  smaller  variance.  In  the  field  tests,  it 
was  possible  to  directly  verify  that  the  voltage  output  from  the  Macrodyne 
frequency  counter  was  proportional  to  the  Doppler  frequency  observed  on  an 
oscilloscope  trace.  We  thus  have  confidence  that  the  Macrodyne  output 
accurately  reflects  the  current  velocity.  It  was  not  possible  to  verify  the 
TSI  frequency  tracker  output  in  the  field  due  to  signal  drift  and  noise. 

While  the  frequency  tracker  output  occasionally  appears  to  follow  the 
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frequency  counter  output,  it  generally  does  not.  This  record  confirms  the 
field  observation  that  the  Doppler  burst  rate  was  not  high  enough  to  allow 
continual,  accurate  frequency  tracking.  Much  of  the  observed  signal  is  a 
result  of  the  tracker  locking  onto  noise  signals.  The  origin  of  the  velocity 
offset  in  the  frequency  tracker  signal  is  probably  associated  with  selection 
of  the  center  frequency  during  the  initial  manual  location  of  the  Doppler 
signal.  When  a  signal  is  not  found,  the  TSI  tracker  tends  to  drift  towards 
this  frequency.  In  this  record,  the  center  frequency  appears  to  correspond  to 
a  velocity  of  about  12  cm/s. 

The  Doppler  burst  rate  or  data  rate  was  obtained  by  counting  the  number  of 
times  the  Macrodyne  output  changed  each  second.  This  rate  is  quite  variable, 
as  shown  in  Figure  10  for  record  199F.  The  mean  data  rate  for  this  record  was 
10  Hz  with  a  range  of  A  to  18  Hz.  The  mean  current  velocity  during  this  time 
was  lA .A  cm/ s ,  which  gives  a  mean  spatial  resolution  of  1.3  cm.  The  spatial 
resolution  varies  due  to  the  random  distribution  of  Doppler  bursts.  Given  a 
record  of  the  data  rate,  it  is  possible  to  select  periods  during  which  the 
data  rate  is  relatively  high  if  an  analysis  is  desired  of  velocities  with  a 
smaller  spatial  resolution. 

The  LDL  laser  diode  unit  was  not  as  well  aligned  as  the  Hitachi  unit.  The 
mean  data  rates  obtained  from  this  unit  are  therefore  considerably  lower.  The 
maximum  data  rate  obtained  from  this  unit  is  just  as  high  as  that  obtained 
from  the  Hitachi.  This  is  because  the  Doppler  burst  rate  should  be  about  the 
same  for  both  units.  The  LDL  rate  is  lower  because  the  signal  strength  was 
lower  and  more  easily  obscured  by  noise.  When  the  noise  level  decreased,  the 
LDL  unit  gave  a  data  rate  just  as  high  as  the  Hitachi. 

Comparison  of  DLDV  and  MRC  Results 

An  extensive  comparison  of  the  velocity  data  measured  simultaneously  with 
the  two  experimental  DLDV  units  and  one  of  the  MRCs  was  presented  by  McPhee 
(1985b).  A  few  examples  are  excerpted  here  from  that  report  to  demonstrate 
the  performance  of  the  DLDV  in  the  Arctic  environment.  Attempts  are  also  made 
to  interpret  and  correlate  the  velocity  data  measured  by  the  two  systems. 

Figure  11  shows  a  comparison  of  the  1-second  averaged  time  series  (lower), 
the  spectra  (upper)  and  the  squared  coherency  (middle)  of  the  horizontal 
velocity  component.  Note  that  the  MRC  time  series  is  shifted  downward  A  cm/ s 
to  facilitate  comparison.  The  DLDV  unit  for  measuring  this  component  was 
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Figure  10.  Data  Rate  at  1-second  Intervals  During  Record  199F 


equipped  with  the  LDL  diode  laser,  and  the  Macrodyne  counter  was  used  to 
process  the  Doppler  signal.  From  the  two  time  series,  a  high  degree  of  simi¬ 
larity  between  the  traces  is  noted,  especially  for  the  medium-  to  low-frequency 
components.  As  shown  by  McPhee  (1985b),  the  horizontal  velocity  components  of 
all  the  runs  generally  compare  quite  well.  On  the  average,  the  mean  and 
variance  are  about  2.4%  higher  and  5%  lower  for  the  DLDV  data  than  for  the  MRC 
data  (Tables  2  through  5,  McPhee,  1985b).  This  is  consistent  with  the  fact 
that  the  MRC  was  50  cm  closer  than  the  DLDV  to  the  bottom  of  the  ice  floe 
within  the  oceanic  boundary  layer. 

The  time  series  of  the  DLDV  data  shows  a  "blocky"  formation,  which  is 
caused  by  a  slight  misalignment  of  the  two  beams  of  the  LDL  unit  as  discussed 
in  Section  5.1.1.  The  misalignment  effectively  reduces  the  data  rate  and, 
thus,  the  spatial  resolution.  As  a  result,  the  dynamic  response  of  this  DLDV 
unit  is  closer  to  that  of  the  MRC.  This  is  demonstrated  by  the  reasonable 
match  of  the  two  spectra.  Apparently,  the  "blocky"  formation  has  little 
effect  on  the  spectral  density  for  the  frequency  range  shown  in  the  figure. 

At  the  high-frequency  end,  both  velocity  spectra  tend  to  follow  approximately 

-5/3 

an  f  '  power  law. 

Figure  12  shows  the  time  series,  spectra  and  squared  coherency  of  the 
velocity  component  at  45®  in  the  vertical  plane  of  the  horizontal  component. 

For  this  configuration, the  DLDV  equipped  with  the  Hitachi  laser  diode  was 
used,  and  the  output  was  fed  to  the  Macrodyne  counter  (channel  0)  and  to  the 
TSI  tracker  (channel  1).  Here,  we  do  not  observe  a  "blocky"  formation,  as  the 
DLDV  with  the  Hitachi  diode  always  stayed  aligned  during  the  experiment.  In 
general,  the  time  series  from  both  channels  0  and  1  are  highly  similar; 
however,  the  magnitude  of  the  velocity  component  measured  with  the  TSI  tracker 
is  consistently  about  20%  lower  than  that  measured  with  the  Macrodyne 
counter.  Part  of  the  reason  is  the  presence  of  a  DC  offset  voltage  that 
introduces  a  small  error  in  the  frequency-voltage  conversion.  It  is  inherent 
to  the  tracker  that  it  would  not  work  well  for  low  burst  rates,  which  is 
confirmed  by  the  field  results. 

Let  us  now  focus  on  the  comparison  of  the  results  measured  with  the 
DLDV/Macrodyne  (Channel  0)  and  the  MRC.  For  the  45®  configuration,  the 
velocity  component  consists  of  both  the  vertical  and  horizontal  components. 
Because  the  MRC  was  about  2  m  and  the  DLDV  2.5  m  below  the  bottom  of  the  ice 
floe,  the  influence  of  the  bottom  topography  would  be  more  pronounced  on  the 
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vertical  velocity  component  measured  with  the  MRC  than  with  the  DLDV.  This  is 
reflected  from  the  noticeably  lower  energy  content  in  the  medium-  to 
low-frequency  components  measured  with  the  DLDV  than  in  those  measured  with 
the  MRC,  which  is  different  from  the  comparison  of  the  horizontal  component 
(Figure  11).  Comparison  of  the  results  of  a  similar  configuration  for  other 
runs  shows  essentially  the  same  trend  (McPhee,  1985b).  For  all  the  runs  with 
the  DLDV  set  at  the  45®  configuration,  the  mean  and  variance  are  10%  higher 
and  8.4%  lower  for  the  DLDV  data  than  for  the  MRC  data,  again  consistent  with 
the  relative  placement  of  the  clusters  in  the  boundary  layer. 

The  above  difference,  resulting  from  the  relative  placement  of  the  DLDV 

and  MRC  in  the  boundary  layer,  may  also  be  identified  from  the  spectra  of  the 

45°  velocity  component.  As  shown  in  Figure  12,  the  spectrum  estimated  from 
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the  MRC  data  follows  an  approximate  f  power  law,  whereas  that  estimated 
from  the  DLDV  data  has  a  more  gentle  slope  of  about  f  A  deficit  in  the 
spectral  density  at  medium-  to  low-frequency  components  effectively  reduces 
the  slope  of  the  spectrum  estimated  from  the  DLDV  data.  A  similar  trend  is 
consistently  observed  for  the  velocity  spectra  for  all  the  45°  configurations 
(McPhee,  1985b). 

The  method  reported  by  Bell  (1983),  as  described  in  Section  4.2,  was  then 
used  to  estimate  the  velocity  spectrum  of  the  DLDV  data.  In  performing  the 
spectral  analysis  of  the  DLDV  data,  only  the  data  points  that  are  updated 
(having  values  differ  from  those  sampled  at  the  previous  time  step)  are  used. 
Such  a  scheme  results  in  a  time  series  sampled  at  irregular  time  intervals 
that  are  multiples  of  the  constant  sampling  time  interval,  which  may  lead  to 
reduction  in  the  spectral  aliasing.  As  demonstrated  by  Masry  et  al.  (1978), 
aliasing  is  avoided  when  the  signal  is  sampled  at  random  time  intervals  . 

Figure  13  shows  a  spectrum  estimated  by  the  method  of  Bell  (1983),  with 

the  cutoff  (Nyquist)  frequency  extended  to  12  Hz,  corresponding  to  the  time 

series  of  channel  0,  record  199G.  As  can  be  seen  in  the  figure,  the  spectrum 

displays  an  f  ^  power  law  extending  more  than  one  decade  from  about  0.03  to 
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1  Hz .  Beyond  1  Hz,  the  slope  of  the  spectrum  reduces  to  about  f  .  The 
method  of  Bell  (1983)  apparently  reduces  but  does  not  remove  entirely  the 
aliasing,  as  the  time  series  was  sampled  at  irregular  but  not  random  time 
intervals.  As  a  result,  a  slight  flattening  of  the  slope  remains  at  the 
high-frequency  end.  Comparison  of  the  spectra  of  the  DLDV  data  shown  in 
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Figure  13.  Velocity  Spectra  Estimated  from  the  DLDV  Data, 
Run  199G-ChO,  Using  the  Method  of  Bell  (1983) 
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Figures  12  and  13  shows  that  both  display  an  f  ^  regime  at  the  high-frequency 
end.  Additional  work  is  required  to  resolve  the  question  of  the  most  appropri¬ 
ate  prefiltering  required  to  obtain  representative  spectra  from  the  DLDV  signal. 

5.1.3  DLDV  Performance  in  MIZ  Environment 

Two  experimental  DLDV  systems  were  successfully  deployed  during  MIZEX  84. 
Both  systems  were  operative  throughout  the  three-week  deployment;  the  one 
equipped  with  the  Hitachi  laser  diode  outperformed  the  one  equipped  with  the 
LDL  laser  diode.  A  total  of  several  hours  of  high-resolution  turbulent  flow 
data  was  recorded  under  an  ice  floe.  The  DLDVs  were  calibrated  against  a  set 
of  ducted  current  meters  placed  in  close  proximity,  and  good  agreement  of  the 
results  was  observed  (see  Section  5.1.2).  Our  assessment  of  the  performance 
of  the  DLDV  in  the  MIZ  environment  is  given  in  this  subsection. 

The  observed  data  rates  were  compatible  with  those  observed  in  the  open 
ocean  with  low  particle  densities  [l] .  However,  the  mean  current  in  the  MIZ 
oceanic  boundary  layer  under  an  ice  floe  was  usually  below  10  cm/s.  On  the 
average,  for  a  sampling  rate  of  48  H?,  a  spatial  resolution  of  between  1  to 
2  cm  for  the  DLDV  data  was  obtained  during  these  tests,  representing  an 
improvement  in  resolution  of  at  least  one  order  of  magnitude  over  mechanical 
turbulence  measurement  techniques. 

The  feasibility  of  the  development  of  a  portable  DLDV  for  deployment  in 
the  Arctic  environment  was  successfully  demonstrated  during  MIZEX  84.  Our 
experience  suggested  several  ways  to  improve  the  performance  characteristics 
of  the  DLDV  for  applications  in  the  Arctic  environment,  as  listed  below.  A 
few  of  these  improvements  were  incorporated  into  the  design  of  the  prototype 
units  subsequently  deployed  during  AIWEX  85. 

(1)  The  laser  power  increases  with  decreasing  temperature.  It  is  essential 
to  use  a  laser  diode  equipped  with  a  photodiode  monitor  to  maintain 
constant  laser  power  regardless  of  temperature  (e.g.,  Hitachi  HL7801 
series)  . 

(2)  Power  levels  used  during  these  tests  were  conservative  to  minimize  the 
chances  of  damage  to  the  lasers.  A  laser  power  of  5  mW  appears  to  be 
adequate,  but  higher  power  would  enhance  the  scattered  light  intensity 
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and  thus  increase  the  burst  rate  or  allow  a  larger  separation  between 
the  tip  of  the  fiberglass  rod  and  the  focal  volume  to  minimize  flow 
disturbance . 

(3)  The  optical  components  should  be  treated  with  anti-reflective  coatings 
to  increase  laser  light  transmission  through  the  transmission  optics 
and  to  reduce  light  reflection  from  optical  surfaces.  It  is  preferable 
to  use  a  single-element  beam  splitter  to  maintain  alignment  of  the 
optics  induced  by  large  temperature  variations. 

(4)  The  processors  were  designed  for  a  relatively  high  burst  rate  and 
high-speed  flows.  Due  to  the  low  burst  rate,  it  is  necessary  to  use  a 
counter  rather  than  a  frequency  tracker  for  processing  the  Doppler 
signals.  Improved  filtering  and  Doppler  burst  discrimination  could 
increase  the  data  rate  and,  hence,  improve  the  spatial  resolution  of 
the  instrument  by  a  factor  of  2  or  better. 

5.2  AIWEX  85 

Three  prototype  HRCs,  each  consisting  of  a  three-axis  DLDV  system  and  a 
high-frequency  temperature/conductivity  probe  pair,  were  deployed  during 
AIWEX  85.  The  objectives  were  to  test  the  system  under  field  conditions  and 
assess  the  system's  reliability,  robustness  and  performance  in  the  Arctic  and 
to  measure  the  high-frequency  turbulent  fluxes  of  momentum,  heat  and  salinity 
in  the  oceanic  boundary  layer  under  ice. 

5.2.1  AIWEX  85  Site  and  Setup 

Observations  of  the  turbulent  boundary  layer  beneath  sea  ice  were  made 
with  the  DLDV  system  during  AIWEX  85,  which  was  carried  out  during  March  and 
April  at  a  camp  in  the  Beaufort  Sea  approximately  200  miles  northwest  of 
Prudhoe  Bay,  Alaska.  The  DLDV  clusters  were  deployed  through  first-year  sea 
ice  with  a  thickness  of  165  cm.  The  deployment  site  was  well  over  100  m  from 
any  pressure  ridge  that  could  perturb  the  boundary  layer  structure.  The  DLDV 
deployment  occurred  during  the  period  from  March  26  to  April  25.  Air  tempera¬ 
ture  during  this  time  ranged  from  -30  to  -17®C,  while  wind  speeds  were  less 
than  5  m/s.  The  boundary  layer  under  these  conditions  is  well-mixed  with  a 
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layer  of  essentially  constant-salinity,  constant-temperature  sea  water  in 
equilibrium,  with  the  sea  ice  extending  down  to  a  depth  of  20  to  30  meters. 

The  DLDV  clusters  were  mounted  on  a  rigid  frame  and  lowered  through  a 
100-cm-diameter  hole  drilled  through  the  ice.  The  orientation  of  the  frame 
was  observed  relative  to  local  magnetic  north. 

Each  DLDV  cluster  consists  of  three  DLDV  transmitter/receiver  pairs.  These 
are  oriented  to  observe  three  orthogonal  components  of  the  flow  field  and  thus 
obtain  a  complete  description  of  the  turbulence  statistics.  In  addition,  each 
three-axis  cluster  has  a  high-resolution  temperature  and  conductivity  probe 
mounted  within  2  cm  of  the  DLDV  observation  volume.  Two  clusters  were  de¬ 
ployed,  one  at  2  and  one  at  4  m  beneath  the  lower  surface  of  the  ice.* 

Each  cluster  included  a  temperature/conductivity  probe  pair.  These  probes 
were  calibrated  at  the  beginning  and  end  of  the  experiment.  Calibration  of 
the  conductivity  probe  was  accomplished  with  a  Beckman  conductivity  measure¬ 
ment  cell.  The  conductivity  of  a  sample  of  snowmelt  and  sea  water  measured  on 
the  Beckman  instrument  was  used  to  calibrate  the  conductivity  probe.  The 
temperature  calibration  was  carried  out  in  two  buckets  of  snow  slush  formed 
with  fresh  water  and  with  sea  water.  The  temperature  of  these  baths  was 
observed  with  a  thermometer  accurate  to  within  0.01®C.  Calibration  was  done 
at  a  gain  of  1,  and  the  measurements  were  made  at  a  gain  of  100  (with  the  DC 
voltage  offset).  Note  that  the  temperature/conductivity  probes  were  designed 
to  measure  differentials,  hence  the  fluxes  but  not  the  absolute  values  of  sea 
water  temperature  and  conductivity,  and  calibration  was  not  conducted  in  a 
well-controlled  environment  in  the  field.  Therefore,  the  absolute  temperature 
and  conductivity  cannot  be  determined  accurately.  We  may,  however,  infer  the 
absolute  temperature  from  the  Sea-Bird  temperature  gauge  deployed  adjacent  to 
the  DLDV-T-C  clusters.  Since  both  the  voltage  outputs  of  the  gauges  are 
linearly  proportional  to  the  temperature  and  conductivity  of  the  sea  water 
within  the  measurement  range  of  interest,  the  relative  change  of  these  physi¬ 
cal  quantities  may  be  accurately  measured  at  higher  gains  (with  the  sensi¬ 
tivity  of  the  gauges  determined  at  the  gain  of  1).  At  a  gain  of  100  and  a 
12-bit  resolution  (10  V  full  scale),  the  gain  may  resolve  minute  changes  of 
0.002“C  and  2  Jjmho/cm. 

*0ne  of  the  processors  was  damaged  due  to  a  power  surge  at  camp. 


The  mean  flow  direction  beneath  the  ice  was  observed  using  a  drifter  such 
as  a  piece  of  gum  or  Tye-wrap.  The  two  clusters  were  then  oriented  so  that 
the  cluster  faced  into  the  current  at  an  angle  of  30  to  45“  to  the  left  of  the 
current.  Both  clusters  were  oriented  in  the  same  direction.  An  observation 
was  then  made  of  mean  current,  temperature  and  conductivity  for  a  period  of 
60  seconds.  The  mean  values  of  temperature  and  conductivity  voltages  were 
input  as  offset  voltages.  The  variations  from  these  reference  values  were 
then  multiplied  by  a  gain  factor  of  100  for  output  to  the  data  file.  The  mean 
update  was  also  used  to  orient  the  clusters,  so  that  the  velocity  and  the 
number  of  data  points  obtained  from  each  cluster  were  similar,  and  to  estimate 
the  data  rate  to  be  expected  . 

The  DLDV  clusters  were  checked  for  alignment  in  air  before  deployment. 

Any  variations  from  alignment  were  corrected,  and  the  system  was  deployed 
just  beneath  the  surface  of  the  water.  A  squeeze  bottle  was  then  used  to 
remove  any  air  bubbles  collecting  around  the  lenses  and  receiver  collimation 
tubes.  The  system  was  then  checked  for  Doppler  signals  using  an  oscillo¬ 
scope.  Preliminary  tests  indicated  that  it  was  necessary  to  readjust  the  gain 
settings  on  each  filter  band  of  the  system  to  match  the  processor  amplifier  to 
each  DLDV  pair.  It  was  necessary  to  disable  the  two  highest  band-pass  filters 
on  the  processor  units  to  remove  a  high-amplitude  radio  signal  from  a  radio 
beacon  used  in  the  camp.  The  lowest  filter  band  was  also  disabled,  as  the 
processor  was  unable  to  filter  out  noise  signals  at  frequencies  in  the  range 
from  256  to  1000  Hz.  Note  that  disabling  the  lowest  frequency  band  increases 
the  threshold  speed  to  about  0.34  cm/s  based  on  the  calibration  of 
2960  Hz/(cm/s) . 

The  sea  water  encountered  was  very  clear;  objects  were  visible  to  depths 
of  20  meters.  This  was  reflected  in  a  fairly  low  Doppler  burst  rate.  While 
it  was  not  possible  to  count  directly  the  Doppler  burst  rate,  it  appeared  to 
be  in  the  range  of  2  to  5  Hz.  This  is  very  similar  to  burst  rates  observed 
during  MIZEX  in  the  Greenland  Sea  in  July  of  1984.  Observations  were  made  at 
a  variety  of  burst  rates  ranging  from  2  Hz  to  25  Hz.  The  spectra  of  turbulent 
motions  observed  by  sampling  at  these  two  rates  overlap,  so  the  sampling  rate 
was  alternated  from  one  sample  to  the  next.  First,  a  sample  would  be  taken  at 
25  Hz  for  a  period  of  300  seconds,  then  a  sample  at  2  Hz  for  3000  seconds.  In 
this  way  the  entire  spectrum  of  turbulent  motions  was  observed. 


48 


TR-352/11-85 


Data  files  were  stored  in  300-kilobyte  files  on  a  RAM  disk.  These  files 
were  transferred  into  1  megabyte  files  on  the  COMPAQ  plus  computer  hard  disk 
memory  using  the  MS-DOS  pipeline  utility.  This  process  takes  about  30  minute 
for  a  typical  file.  Primary  processing  of  the  raw  data  was  then  performed 
using  a  program  called  CULL  (See  Section  4.2),  and  the  results  were  stored  on 
floppy  disk. 

5.2.2  Results  and  Discussion 

A  total  of  28  records  of  velocity,  temperature  and  conductivity  was 
obtained  at  depths  of  2  and  4  m  beneath  the  bottom  of  the  ice.  The  first 
15  records  were  processed  using  faulty  assumptions  about  the  way  in  which  the 
binary  data  file  was  stored  and  are  thus  of  only  minimal  value.  In  addition, 
the  filters  on  the  processing  system  were  not  properly  disabled  until  record 
17.  Table  4  lists  the  last  ten  records,  19  through  28,  along  with  the  sample 
rate,  data  rate  and  mean  velocity.  The  AIWEX  85  data  have  been  partially 
analyzed.  Some  of  the  results  are  presented  and  discussed  in  this  subsection 


Table  4. 

AIWEX 

85  Data  Summary 

Depth  Rate  Period 

Data 

Rates:  U1 ,  U2,  U3 

Speed 

Heading 

Record  Date  (m)  (Hz)  (s) 

(Hz) 

(cm/ s) 

(degrees) 

5 

1200 

2 

3000 

5 

240 

2 

3000 

2  3nj0 


2  3000 


Time  Series 

The  record  length  of  the  DLDV  runs  ranged  from  240  to  3000  seconds.  In 
only  one  of  the  observations  in  the  boundary  layer  did  the  variation  in 
temperature  and  conductivity  exceed  0.01®C  and  10  vmho/cm,  indicating  a 
well-mixed  boundary  layer  for  a  time  scale  of  up  to  one  hour  during  AIWEX  85 
(March  and  April  1985).  Measurements  by  other  AIMEX  participants  (McPhee , 
private  communication)  show  essentially  constant  values  of  these  quantities 
over  periods  comparable  to  that  of  the  DLDV  runs.  However,  they  show 
noticeable  long-term  variations  in  these  quantities.  Figures  14  through  16 
show  the  typical  time  series  of  temperature  and  conductivity  variations  and 
the  U1  velocity  component  at  the  2-m  depth.  The  temperature  time  series  is 
practically  uneventful,  whereas  the  conductivity  time  series  shows  only 
high-frequency  fluctuations.  Only  the  U1  velocity  component  shows  noticeable 
temporal  variations  in  both  high-  and  low-frequency  ranges. 

Figures  17  through  19  are  another  set  of  time  series  of  the  same 
quantities.  Figure  18  shows  that  the  temperature  time  series  at  the  4-m  depth 
displays  a  wave-like  feature  with  a  period  of  about  10  minutes  and  a  wave 
height  of  0.03'’C.  This  sharp-crested  wave-like  feature  of  low  amplitude  is 
the  only  one  observed  among  all  the  time  series.  A  corresponding  feature  is 
not  observed  in  the  conductivity  and  velocity  time  series  from  the  same  run. 
Again,  only  high-frequency  fluctuations  are  observed  in  the  conductivity  time 
series,  whereas  the  velocity  components  show  a  relatively  wide  range  of 
temporal  variation. 

Mean  Current  Speed  and  Heading 

The  mean  current  speed  for  the  10  records  listed  in  Table  4  ranges  from 

about  4  to  12  cm/ s .  The  mean  current  speed  at  the  4-m  depth  is  consistently 

higher  than  that  at  the  2-ra  depth,  indicating  that  the  DLDVs  are  within  the 

boundary  layer  in  which  a  vertical  shear  is  present.  For  all  ten  records,  the 

mean  current  shear,  defined  as  (U  ,  -  U  „  )/U  ,  in  a  dimensionless 

z=4m  z=2m  z=4m 

form  where  ^^=2^  ^z=4m  mean  current  speeds  at  the  2-  and  4-m 

depths,  is  0.07.  In  addition,  the  mean  headings  of  the  current  at  the  two 
depths  also  differ  by  about  14®  with  the  current  vector  rotating  in  the  clock¬ 
wise  direction  as  the  depth  increases.  This  turning  of  the  current  vector 
corresponds  to  the  spiral  structure  in  the  Ekman  layer  (Ekman,  1905).  The 
effect  of  the  Ekman  spiral  was  also  observed  by  Hunkins  (1975)  during  the 
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Arctic  Ice  Pynamics  Joint  Experiment  pilot  experiment.  The  angle  of  deflec¬ 
tion  between  the  two  current  vectors  at  the  2-  and  4-m  depths  is  generally 
larger  for  our  measurements  than  those  measured  by  Hunkins  (1975).  Further¬ 
more,  our  results  show  that  the  angle  of  deflection  tends  to  increase  as  the 
current  speed  decreases.  For  example,  the  angle  of  deflection  for  the  run 
with  the  highest  current  speed  (record  24,  U  ,  =  12.3  cm/s)  differs  by  12°, 

whereas  that  for  the  lowest  speed  (record  25,  U  ,  =  4.7  cm/s)  differs  by 

z— 4m 

28°.  The  same  trend  holds  for  the  runs  with  the  next  highest  and  lowest 
current  speeds . 

Intuitively,  the  bottom  topography  of  the  ice  floe,  such  as  the  presence 
of  pressure  keels,  would  have  a  relatively  strong  local  influence  on  the 
headings  of  the  current  at  these  two  depths. 

Spatial  Resolution 

From  the  results  given  in  Table  4,  we  estimate  the  spatial  resolution, 
which  is  equal  to  the  mean  current  speed  divided  by  the  mean  data  rate  of  the 
U1  velocity  component  (Figure  20).  Also  plotted  in  the  figure  is  the  spatial 
resolution  estimated  from  the  MIZEX  84  data.  Figure  20  clearly  demonstrates 
the  trend  of  improving  spatial  resolution  with  increasing  sampling  rate.  In 
an  Arctic  environment,  a  spatial  resolution  of  2  cm  may  be  achieved  for  a 
sampling  rate  of  25  Hz  and  higher.  Therefore,  based  on  the  above  findings,  at 
higher  particle  densities  in  water  bodies  other  than  the  Arctic  Ocean  together 
with  an  increased  sampling  rate,  a  spatial  resolution  of  the  DLDV  down  to  1  mm 
is  achievable  . 

The  spatial  resolution  estimated  from  the  AIWEX  85  data  could  be  further 
reduced  partly  by  optimizing  the  design  of  the  counter  electronics  of  the 
three-channel  FLOW  counter  processor.  As  a  trade-off  to  the  multi-channel 
design,  each  of  the  three  channels  of  the  processor  makes  only  one  5/8 
comparison  during  each  sampling  cycle  (See  Section  2.3.1).  After  detecting  a 
burst  signal  that  does  not  pass  the  validation  criterion,  the  processor  would 
not  look  for  another  burst  signal  even  if  there  is  enough  time  to  perform  such 
comparison  during  the  same  cycle.  For  the  single-channel  Macrodyne  processor 
used  in  MIZEX  84,  however,  the  5/8  comparison  was  made  at  a  higher  rate  than 
the  sampling  rate.  Therefore,  a  higher  data  rate  or  finer  spatial  resolution 
was  attained  with  the  Macrodyne  processor  than  with  the  FLOW  processor,  as  the 
spatial  resolution  is  inversely  proportional  to  the  data  rate.  Software  and 
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Figure  20.  Estimated  Spatia 


hardware  modifications  may  be  incorporated  into  the  FLOW  processor  to  perform 
an  essentially  continuous  5/8  comparison  to  increase  the  data  rate  and,  thus, 
to  reduce  the  spatial  resolution. 

Wavenumber  Spectra 

The  wavenumber  spectra  of  the  velocity  components  were  estimated  using  the 

method  of  Bell  (1983).  Figure  21  is  a  composite  wavenumber  spectrum  of  the  U1 

component  measured  at  2-  and  A-m  depths  for  Run  No.  AIWEX  24.  The  abscissa 

and  ordinate  are  the  wavenumber  k  in  cm  ^  and  the  normalized  velocity 

spectrum.  Because  the  individual  spectra  are  essentially  the  same,  the 

composite  spectrum  is  presented  to  improve  the  accuracy  of  the  spectrum.  The 

corresponding  wavenumber  spectrum  for  the  U3  component  is  shown  in  Figure  22 . 

-5/3 

Both  spectra  tend  to  show  a  k  trend  for  almost  one  decade  at  wavenumbers 
beyond  the  first  spectral  peak.  The  spectral  density  at  higher  wavenumbers 
has  a  more  gentle  slope  and  eventually  levels  off. 

The  leveling  off  of  the  spectrum  is  mainly  due  to  the  presence  of  infre¬ 
quent  spikes,  which  contaminate  the  signals,  and  to  limitations  in  preprocess 
filtering,  which  tend  to  enhance  high-frequency  spectral  power.  As  described 
in  Section  5.2.1,  one  of  the  sources  of  the  spikes  was  a  high-amplitude  radio 
signal  from  a  radio  beacon  used  in  the  camp.  It  was  discovered  that  the 
three-channel  FLOW  counter  processor  that  uses  the  filter  band  concept 
(Section  2.3.1)  had  not  been  optimized  for  noise  isolation  or  filtering. 
Band-pass  filtering  of  the  signals  to  remove  noise  outside  the  frequency  bands 
of  interest  would  definitely  improve  the  performance  of  the  FLOW  processor. 

In  reducing  the  raw  data,  the  large  spikes  were  removes  through  a  persistence 
check  described  in  Section  5.2.1.  Only  those  spikes  that  have  amplitudes  less 
than  2  to  3  times  the  variance  of  the  velocity  were  not  removed.  Although  the 
energy  content  of  these  spikes  is  insignificant  compared  with  the  total  kinetic 
energy,  the  spikes  do  mask  the  low-level  spectral  components  at  the  high- 
wavenumber  end,  especially  when  the  flow  is  not  strongly  turbulent  as  was 
encountered  during  AIWEX  85.  A  better  spike  removal  algorithm  will  be  required 
to  discriminate  the  relatively  small  spikes  from  the  signals.  For  future 
field  experiments,  however,  emphasis  should  be  made  to  improve  the  capability 
of  the  processor  in  noise  isolation  and  suppression.  Several  recommendations 
are  given  in  Section  7  for  further  system  improvements  of  the  DLDV  sensors  and 
processor  for  field  deployment  under  hostile  environmental  conditions. 
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During  the  first  12  months  of  the  Phase  TI  research  and  development  effort, 
we  participated  in  two  Arctic  field  experiments,  using  two  experimental  models 
of  DLDVs  during  MIZKX  84  and  three  prototype  HRCs  of  velocity-temperature- 
conductivity  (V-T-C)  sensors  during  AIWEX  85  for  turbulence  measurements  under 
ice  floes.  The  experimental  DLDV  units  and  the  sensors  and  processors  of  the 
prototype  HRCs  were  successfully  developed,  fabricated,  tested  and  calibrated 
prior  to  deployment.  The  results  from  MIZEX  84  were  also  used  to  revise  the 
design  of  the  prototype  HRC  units  before  their  deployment  during  AIWEX  85. 

In  essence,  we  have  demonstrated  even  in  an  Arctic  environment  that  the 
DLDV  is  a  viable  oceanographic  instrument  for  high-frequency  and  high- 
resolution  turbulence  velocity  measurements.  The  field-proven  DLDV  transmitter 
and  receiver  have  high  potential  to  be  developed  into  a  commercial  product 
without  major  modifications.  At  present,  we  have  partially  analyzed  the  AIWEX 
85  data.  Based  on  our  experience  gained  from  the  field  deployment  and  the 
analysis  performed  so  far,  we  summarize  several  of  our  important  findings: 

(1)  The  field  worthiness  of  the  prototype  DLDV  unit,  including  the 

transmitter  and  receiver,  has  been  proven  in  terms  of  its  portability, 
robustness  and  performance.  The  DLDV  units  were  not  only  trouble-free 
during  the  field  deployment,  but  they  also  performed  satisfactorily  as 
projected  and  survived  the  shipping/handling  and  the  deployment  in  the 
hostile  Arctic  environment.  The  DLDV  unit  in  its  current  configuration 
(extendable  to  20  mW  of  laser  power)  may  readily  replace  a  gas  laser 
of  the  same  power  rating  and  the  associated  optics  and  photodetector 
used  in  commercial  LDV  processors  (counters  or  frequency  trackers)  for 
laboratory  and  field  measurements,  thus  achieving  one  of  the  main 
goals  set  for  the  project.  The  main  advantages  of  the  DLDV  over  a 
conventional  LDV  using  a  gas  laser  and  an  assembly  of  optical  and 
electronic  components  are: 

o  Compactness,  portability  and  robustness  for  laboratory  and  field 
deployment . 

o  Integrated  transmitting  and  receiving  modules  to  simplify  the  test 
procedure . 
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o  Low  power  requirement  for  field  deployment  in  remote  locations . 

Six  months  after  AIWEX  85,  we  turned  on  two  transmitters  selected  at 
random  to  learn  that  they  were  still  in  good  working  order  with  no 
realignment  required.  According  to  the  manufacturer,  Hitachi 
Corporation,  the  lifetime  of  the  laser  diode  is  10^  hours  at  25°C . 

One  of  the  units  is  now  under  life  testing  to  determine  its  longevity 
under  normal  operating  conditions. 

(2)  Simultaneous  measurements  of  the  current  velocity  components  using  the 
ducted  current  meters  on  the  medium-resolution  cluster  (MRC)  and  the 
DLDV  units  (McPhee,  1985b)  placed  in  close  proximity  (50  cm  vertical 
separation)  were  conducted  during  MIZEX  84.  Comparison  of  the  results 
shows  generally  good  agreement  between  the  time  series  and  frequency 
spectra,  especially  for  the  horizontal  velocity  component.  For  the 
45®  velocity  component,  there  is  a  deficit  in  the  spectral  density  of 
the  DLDV  spectrum  at  medium  to  low  frequencies.  The  deficit  is 
apparently  due  to  the  more  pronounced  influence  of  the  ice  bottom 
topography  on  the  vertical  velocity  component  measured  with  the  MRC 
than  with  the  DLDV,  as  the  MRC  is  0.5  m  closer  to  the  ice  bottom  than 
the  DLDV.  The  capability  of  measuring  high-frequency  velocity 
components  (f  >0.5  Hz)  using  the  DLDV  is  demonstrated. 

(3)  From  the  results  of  MIZEX  84  and  the  experience  gained  during  that 
field  experiment,  we  have  successfully  demonstrated  that  the  DLDV  is  a 
viable  instrument  deployable  in  the  Arctic  environment.  Our 
assessment  of  its  performance  during  MIZEX  84  indicates  several 
achievable  improvements  (see  Section  5.1.3)  that  were  subsequently 
incorporated  into  the  design  of  the  prototype  units  (Section  2) 
deployed  during  AIWEX  85. 

(4)  The  time  series  measured  during  AIWEX  85  at  2-  and  4-m  depths  below 
the  bottom  of  an  ice  floe  shows  that  the  variations  in  temperature  and 
conductivity  generally  do  not  exceed  0.01®C  (except  in  one  case)  and 
10  vmho/cm  for  a  time  scale  up  to  an  hour  or  so,  indicating  a 
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well-mixed  boundary  layer  during  April  and  May  of  1985.  This  is 
consistent  with  measurements  by  other  participants  over  time  periods 
comparable  to  that  of  the  DLDV  runs,  although  noticeable  long-term 
variations  were  reported. 

(5)  In  both  locations  where  the  two  field  experiments  were  conducted, 
there  were  adequate  scatterers  in  the  Arctic  Ocean  water  to  achieve 
the  target  spatial  resolution  of  2  cm  or  better  for  velocity  measure¬ 
ment.  Our  results  show  that  the  spatial  resolution  improves  with  the 
data  rate.  A  2-cm  spatial  resolution  was  achieved  with  a  sampling 
rate  of  about  25  Hz  or  higher.  In  other  oceans  and  water  bodies  where 
natural  scatterers  are  plentiful,  a  spatial  resolution  of  the  order  of 
1  mm  is  achievable  when  off-axis  detection  is  employed  at  similar  or 
higher  sampling  rates. 

(6)  The  wavenumber  spectra  of  the  velocity  components  estimated  from  the 

-5/3 

AIWEX  85  data  tend  to  show  the  presence  of  an  f  regime.  However, 
the  spectral  densities  at  high  frequencies  were  contaminated  by  the 
presence  of  infrequent  noise  spikes  induced  by  a  radio  beacon.  As  a 
result,  the  spectrum  at  the  high-frequency  end  levels  off,  as  opposed 
to  a  slight  flattening  of  the  slope  of  the  spectrum  of  the  MIZEX  84 
data.  Optimization  of  the  FLOW  counter  processor  for  better  noise 
isolation/filtering  would  be  required  to  reject  most  of  the  noise 
spikes  . 

(7)  As  a  potential  oceanographic  instrument  for  velocity  measurements,  the 
DLDV  has  several  advantages  over  mechanical  sensors: 

o  High-frequency  response  limited  only  by  the  number  density  of 
scatterers  . 

o  Small  sampling  volume  and  high  spatial  resolution  (down  to 
millimeters) . 

o  Linear  and  absolute  calibration  (no  drift). 
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o  Very  low  threshold  speed  (<  0.3  cm/s). 


o  Minimal  disturbance  to  the  flow  field  (no  physical  probe  at  the 
sampling  volume) . 

7 .  RECOMMENDATIONS 

Based  on  the  results  of  this  two-phase  investigation,  we  have  concluded 
that  the  DLDV  has  a  high  market  potential  to  become  a  commercial  velocity 
sensor  for  field,  laboratory  and  engineering  applications.  To  further  enhance 
its  performance  for  a  wide  range  of  applications,  we  recommend  several  improve¬ 
ments  to  be  incorporated  into  the  design  of  the  DLDV  and  the  processor.  In 
addition,  suggestions  for  software  development  for  post  analysis  of  the  DLDV 
data  obtained  in  MTZEX  84  and  AlWEX  85  are  given  to  enhance  the  SNR  so  that  we 
may  obtain  accurate  estimates  of  the  high-order  statistics  from  the  data  whose 
high-frequency  components  are  contaminated  with  infrequent  noise  spikes. 

(1)  For  the  current  3-axis  DLDV  system,  we  have  avoided  the  directional 
ambiguity  inherent  to  laser  velocimetry  by  orienting  the  three  DLDV 
units  in  a  special  geometrical  configuration  (Section  3.1).  This 
avoidance  scheme  only  works  in  a  flow  field  with  a  prevailing 
current.  For  general  applications,  it  is  preferable  to  have  a  truly 
directional  sensing  capability.  Several  techniques  are  available  for 
removal  of  directional  ambiguity,  but  none  is  compatible  with  the 
compactness  and  portability  of  the  DLDV  for  field  deployment  (Durst 
and  Zare ,  1974).  During  Phase  II,  we  looked  into  several  innovative 
alternatives.  Among  these,  the  technique  of  imaging  the  focal  volume 
onto  a  position-sensitive  or  lateral-effect  photodiode  seemed  most 
promising.  The  lateral-effect  photodiode  has  two  back  contacts  at 
opposite  ends  of  a  single-axis  model.  The  current  output  at  each  of 
the  contacts  increases  with  the  proximity  of  a  light  spot  imaged  onto 
the  active  strip  of  the  diode.  By  imaging  the  focal  volume  of  the 
DLDV  onto  the  diode,  with  the  plane  of  the  beam  intersection  aligned 
with  the  axis  of  the  diode,  we  may  discriminate  whether  a  particle  is 
entering  from  one  end  or  the  opposite  end  of  the  focal  volume,  as 
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detected  from  the  differential  output  of  the  two  contacts  of  the 
diode,  which  senses  the  direction  of  the  light  scattered  from  the 
moving  particle.  From  the  sum  of  the  two  outputs,  we  reconstruct  the 
Doppler  signals  for  sensing  the  speed  of  the  particle.  Currently,  we 
are  testing  this  technique  and  assessing  the  feasibility  of  its 
implementation  to  the  current  DLDV  design  under  the  sponsorship  of 
DOE.  If  successful,  we  will  have  the  option  to  increase  the  laser 
power  and  to  operate  the  DLDV  in  a  back-scattering  mode  by  combining 
two  separate  modules  into  a  single  one,  greatly  simplifying  the 
deployment  of  the  instrument. 

(2)  As  described  in  Section  5.2.2,  the  counter  logic  in  the  three-channel 
FLOW  processor  does  not  have  the  provision  to  look  for  more  then  one 
Doppler  burst,  whether  it  is  validated  or  not,  during  each  sampling 
cycle.  This  effectively  reduces  the  data  rate.  This  provision  should 
be  incorporated  through  hardware/software  design  so  that  the  processor 
continuously  searches  for  and  validates  Doppler  bursts  detected  by  the 
photosensor  until  a  validated  burst  is  encountered  and  updated  during 
each  sampling  cycle. 

(3)  During  AIWEX  85,  we  experienced  interference  of  the  DLDV  signals  from 
a  radio  beacon.  As  a  result,  we  had  to  manually  disable  the  highest 
frequency  band  (>  64  kHz)  of  the  DLDV.  The  lowest  band  (256  to 

1000  Hz  corresponding  to  a  current  speed  of  less  than  0.34  cm/s)  was 
also  disabled  to  avoid  noise  contamination  in  this  frequency  range. 
After  this  correction,  we  still  found  infrequent  noise  spikes  that 
contaminated  the  high-frequency  components  of  the  spectra.  To  improve 
the  performance  of  the  DLDV  for  deployment  in  a  "noisy"  environment, 
better  noise  isolation/f iltering  through  careful  shielding  of  the 
cable  and  instrument  will  be  required.  Furthermore,  provisions  for 
switching  off  individual  filter  bands  should  be  implemented  so  that 
the  optimum  selection  of  the  filter  bands  for  noise  rejection  may  be 
made  under  various  flow  conditions. 
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(4)  A  parallel  effort  to  design  a  decontamination  algorithm  to  recognize 
and  remove  the  infrequent  spikes  that  contaminate  the  DLDV  signals 
would  help  to  recover  the  high-frequency  information.  A  similar 
algorithm  has  been  used  successfully  by  Dugan  and  Okawa  (1982)  to 
remove  spikes  contaminated  by  the  conductivity  signals  due  to 
particulates  passing  through  the  sensing  volume. 

(5)  A  part  of  the  results  derived  from  AIWEX  85  and  presented  in  this 
report  has  shown  some  interesting  features  in  the  high-frequency  and 
small-scale  regimes  that  cannot  be  attained  with  conventional  sensors  . 
This  unique  set  of  data  definitely  deserves  careful  examination  in  an 
attempt  to  characterize  the  turbulence  processes  taking  place  in  the 
Arctic  Ocean  boundary  layer  under  ice.  We  recommend  a  detailed  data 
analysis,  including  the  development/adaptation  of  a  spike  decontamin¬ 
ation  algorithm  and  refinement  of  the  statistical  and  spectral 
analysis,  to  derive  statistical  and  spectral  parameters  of  the  HRC 
data  . 
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APPENDIX  A 

DLDV  OPERATING  INSTRUCTIONS 

Software  Version  1.2 

November  8,  1985 

Flow  Industries  Inc. 

Automation,  Instrumentation  &  Controls  Division 


David  R.  Schmidt 


I 


CONNECTIONS  BETWEEN  CLUSTERS  AND  PERSONAL  COMPUTER 


All  control  and  data  signals  are  passed  over  the  GPIB  cable,  which 
connects  the  personal  computer  (PC)  with  the  cluster  units.  The  cables 
themselves  just  connect  all  of  the  connectors  on  all  of  the  devices  in 
parallel,  although  some  practices  should  be  followed: 

1.  The  units  should  be  connected  in  a  linear  fashion. 


If  the  PC  had  been  connected  in  the  middle  of  the  string  of  other 
devices,  it  may  have  worked  most  of  the  time,  but  reflections  on  the 
cable  could  cause  an  increased  error  rate. 

2.  Don't  over  tighten  the  thumbscrews  that  hold  the  GPIB  connectors. 

They  are  only  there  to  keep  the  cable  from  torquing  the  connector 
loose.  If  they  are  tightened  too  much,  the  threaded  sockets  on  the 
chassis  can  come  loose  or  break. 

3.  The  GPIB  cables  are  expensive  and  hard  to  fix:  beware. 

II.  ERROR  RECOVERY 

Due  to  the  nature  of  the  GPIB,  there  are  times  when  all  of  the  cluster 
units  will  need  to  be  reset  to  escape  from  an  error  condition.  Generally, 
when  the  program  is  'hung  up'  in  some  undesired  state,  reset  all  of  the 
cluster  units  and,  if  necessary,  reboot  the  PC.  A  reset  push  button  is 
located  inside  the  front  door  of  the  cluster  unit  on  the  left  side  of  the 
card  cage.  Cycling  the  power  off  for  10  seconds  should  perform  a  reset, 
but  at  times  it  may  be  necessary  to  press  the  reset  button  as  well.  It 
appears  that  some  of  the  purchased  boards  in  the  cluster  units  do  not 
reset  correctly  from  the  reset  button  and  need  to  be  reset  by  turning  off 
the  power . 

III.  RUNNING  THE  DLDV  PROGRAMS 

The  simplest  way  to  run  the  software  that  makes  the  system  work  is  to 
boot  the  PC  from  the  'DLDV  BOOT  DISK.'  This  diskette  will  set  up  the 
device  drivers  for  the  GPIB  interface,  create  a  ram-based  simulated  disk, 
copy  the  latest  version  of  the  operating  program  to  the  ram  disk  and 
perform  a  few  dos  initializations.  There  are  several  files  of  interest: 

1.  c :''DLDV''DLDV.exe  -  the  main  operating  program 

2.  c ••''DLDVNahowdldv .exe  -  inspect  data  and  status  files 

3.  c :N)LDVNevb .exe  -  special  version  of  showdldv 

4.  c :N)LDVM)LDV.sav  -  boot-up  version  of  DLDV  setup 

5.  b:DLDV.sav  -  current  version  of  DLDV  setup 
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DETAILS  OF  EACH  PROGRAM 


A .  DLDV 
1 .  ex 

Exit  from  DLDV.  Saves  status  of  the  program  for  later  loading  with  the  'rp' 
command  or  inspection  with  showdldv  or  sbv. 

2  .  rp 

On  exit  from  the  DLDV  program,  most  of  the  status  information  is  saved  in  the 
file  'DLDV.sav'  .  This  means  that  it  is  possible  to  escape  from  DLDV  and 
restore  your  place  later  by  using  the  'rp'  command. 

3 .  on  cluster_number  (enable  cluster) 

4.  of  cluster_number  (disable  cluster) 

5.  id  cluster  number  (set  GPIB  address  for  cluster) 


Each  cluster  has  two  numbers  associated  with  it.  One  is  the  cluster  number, 
which  determines  the  column  that  a  cluster's  information  is  displayed  in  on 
the  screen.  The  leftmost  column  is  always  cluster  number  1.  The  other  number 
is  set  by  the  hardware  on  the  GPIB  board  in  the  cluster  unit  itself.  This 
allows  a  particular  hardware  unit  to  be  addressed  by  the  PC.  The  'GPIB 
address'  line  on  the  DLDV  status  display  shows  which  cluster  number  maps  to 
which  of  the  three  cluster  boxes.  It  is  recommended  that  the  GPIB  address 
values  be  left  as  1,  2,  and  3  and  that  they  match  the  cluster  numbers.  This 
may  reduce  confusion.  Internally  to  the  DLDV  program,  each  of  the  clusters 
may  be  enabled  or  disabled.  A  disabled  cluster  is  never  read  or  written  to. 

An  attempt  to  read  or  write  to  a  cluster  that  is  either  not  connected  or 
turned  off  will  cause  the  system  to  hang  up. 

6  .  r  t  t  ime 

(set  run  time  for  mu  and  td ) 

7.  fs  freq 

(set  sampling  frequency  for  mu  and  td) 


8.  % 


(set  5/8  ratio  check  error  bandwidth) 


These  commands  set  the  sampling  parameters  for  the  data  collection.  The  rt 
command  is  followed  by  the  runtime  in  seconds.  Fs  command  frequencies  are 
limited  to  2,  4,  3,  10,  20,  and  25  hertz.  Exit  saves  the  runtime  and  sampling 
frequency,  but  not  the  error  bandwidth  specified  by  the  %  command.  In  doing 
the  5/8  ratio  check,  it  is  not  always  desired  to  require  an  exact  5/8  ratio 
and  the  percentage  allowed  is  specified  with  the  %  command.  THE  PERCENTAGE  IS 
NOT  SAVED  ON  EX  NOR  LOADED  WITH  RP . 
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9.  mu  (update  means  for  all  on-line  clusters) 

This  command  performs  two  functions.  It  allows  the 
operator  to  get  a  quick  check  on  the  average  Doppler 
frequencies  from  all  of  the  on-line  clusters  without  having 
to  collect  data  to  disk.  It  also  sets  up  the  temperature 
and  conductance  board  for  normal  operation. 

10.  td  filename  (collect  data  to  disk) 

This  is  the  real  workhorse  function.  It  will  save  the  data 
collected  from  all  on-line  clusters  into  the  file 
specified.  The  file  should  be  on  the  b:  ram  disk  since  the 
delays  caused  by  a  mechanical  disk  are  too  long.  The 
filename  should  be  a  simple  name  with  no  paths  (no 
backslash  characters  allowed)  .  No  screen  updates  are  done 
during  the  collection,  as  the  computer  is  very  busy 
transferring  the  data  and  no  time  is  left  to  run  the 
display.  At  the  end  of  the  sample  period,  the  elapsed  time 
is  displayed  and  the  system  prompt  appears. 

11 .  gt  cluster_number 

12.  gc  cluster  number 

13.  ct  cluster_number 

14.  cc  cluster_number 

Adjust  the  temperature  and  conductivity  calibrations  for 
the  various  clusters.  Gains  are  specified  with  the  single 
digits  0,  1,  2,  or  3,  which  correspond  to  physical  gains  of 
1,  10,  100,  and  1000  respectively. 

15  .  lo  placename 

The  placename  is  up  to  16  characters  long  and  is  saved  with 
the  data  and  status  files.  It  is  not  used  by  the  program, 
so  it  may  contain  any  kind  of  identifying  information 
desired.  The  placename  may  also  be  set  from  the  dos 
command  line.  For  example  the  dos  command 

DLDV  kent  wa  001 

would  run  DLDV  and  set  the  placename  to  'kent  wa  001'. 


16 .  ud 


This  is  useful  for  verifying  that  the  GPIB  connection 
works.  It  merely  transfers  some  data  to  the  on-line  units. 
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17  .  cz  cluster  number 


Clears  the  cluster  specified.  Not  especially  useful  except 
during  software  debugging. 


18  .  end 


Causes  an  immediate  exit  from  the  DLDV  program.  DOES  NOT 
SAVE  STATUS.  Useful  if  the  disk  fills  up  or  you  don't  want 
to  save  the  current  status  for  some  reason. 

V.  MISCELLANEOUS  NOTES 

1.  If  its  on  the  ram  disk  (b:)  and  you  want  to  keep  it,  copy  it  to  the 
hard  disk  (c:)  or  to  a  floppy  (a:)  before  turning  the  machine  off. 

2.  You  won't  get  any  sympathy  for  not  doing  backups  of  both  your  floppies 
and  hard  disk. 

3.  The  system  clock  must  be  set  from  dos . 
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APPENDIX  C 

DESIGN  INFORMATION  FOR  THE  TEMPERATURE /CONDUCTIVITY  CIRCUITRY 

PTC-1 

Programmable 

TEMPERATURE 

CONDUCTANCE 


by 

Paul  Stockton 

Datacon  Electronics 
Albany,  Oregon 
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OCT  30,84 

TEMPERATURE  -  CGMDOCTANCE 
CIRCUIT  DESCRIPTIW: 

The  unit  is  a  single  circuit  board  designed  to  be  physically  and 
electranicadly  compatible  with  an  industry  standard  computer  internal 
buss  structure  known  as  the  Standard  Buss  (STD-BUS) . 

The  prototype  board  performs  three  basic  functions. 

1. )  To  convert  ten^rature  to  a  voltage  acceptable  to 

the  anedog  interface  board. 

2. )  To  convert  the  conductance  of  Scdtwater  to  a  voltage 

acceptable  to  the  analog  interface  board. 

3 .  )  Provide  programmable  control  of  both  temp  &  ccnd 

to  enable  cooputer  control  of  gain  and  offset 
for  both  top  &  cond.  This  capability  provides 
very  high  resolution  over  a  wide  dynamic  range. 

TEMPERATURE;  The  tenp.  circuit  is  a  siaple  resistance  divider 
where  one  resistor  is  the  temperature  sensitive  element  (fast  probe 
model  EP-07)  and  a  precision  load  resistor  (Rl).  Voltage  to  this 
divider  is  provided  by  a  precision  voltage  reference  source.  The 
value  of  Rl  has  been  choosen  to  provide  maxlmim  sensitivity  at  the 
anticipated  tenperature  range  of  -4  degC  to  -flS  degC.  A  very  low 
noise  low  drift  precision  amplifier  is  the  voltage  follower /buffer 
for  the  signal  from  the  thennister.  The  signal  is  then  Low  Pass 
filtered  at  either  25Hz,  50Hz  or  iOOHz  selectable.  The  signal 
is  now  a  voltage  representing  temperature  and  is  controllable 
by  the  computer  through  a  selectable  gain  amplifier  and  offsettable 
by  a  sixteen  bit  digital  to  voltage  converter.  The  possible  gains  are 
XI,  XIO,  XlOO,  and  XIOOO.  Prior  to  the  gain  control,  the  signal  is 
offsettable  from  OVDC  to  '•■lO.OOOOVDC  in  incremants  of  0.0CX)153VDC. 

CX3NDUCTANCE;  Conductance  is  determined  by  RMS  value  of  an  AC 
current  (Isig)  flowing  through  the  sample  solution  being  tested. 

The  conductance  is  equal  to  the  probe  voltage  (Vsig)  divided  by 
the  probe  current  (Isig) .  The  circuit  measures  both  Isig  and 
Vsig,  divides;  Isig/Vsig  to  obtain  a  voltage  proportional  to 
conductance  independent  of  fluctuations  in  Vsig  (within  limits) . 

This  voltage  is  then  processed  through  exactly  the  circuitry  as 
the  temperature  so  control  circuitry  has  matched  freq.  response  and 
phase  shift.  Again,  the  programnable  gain  and  offset  provide  very 
large  dynamic  range  and  resolution. 

NDTEiThe  RMS — >DC  converters  in  the  conductance  circuit,  may  produce 
a  snail  coorelation  error  of  the  temp  vs  cond.  They  also 
limit  freq  response.  These  can  be  eliminated  with  further  development. 
Post  processing  of  the  raw  data  can  remove  this  phase  shift  difference 
If  bbjectionable. 
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TERKENOLOGY  OCT  31,84 


Alp  :  If  gain  of  IP  filters. 

Abp  :  Bandpass  gain  of  the  IP  filters. 

B  :  One  of  two  equations,  "A"  &  "B"  for  filters. 

Cave  :  Averaging  capacitor  for  the  RMS — >DC  converter. 

*Note:  Will  determine  the  freq  res  of  the  RMS — >DC  conv. 
[CIRCUIT  TRADE-OETS:  1.)  Keep  the  freq  of  the  sine  wave 
oscillator  as  hi^  as  possible  to  keep  Cave  as  snail  as 
possible  to  keep  the  response  of  the  conv.  high. 

2 . )  The  freq  is  limited  by  the  ix^t  asps  high  freq  response 
If  the  input  anps  gadn  changes  with  signal  strength, 
non-linearity  will  occur. 

Cext  :  External  capacitor  for  the  IP  filters  UAF14. 

Cfo  :  The  frequency  determining  capacitor  of  the  sine  wave  osc. 

Dout  :  DC  voltage  at  output  of  conductance  divider. 

fo  or  fc  :  -3db  point  at  freq.  fo  of  IP  filters.  (UAF41) 

Freq  or  Fco  :  Frequency  of  Sinewave  Oscillator. 

IRa  :  Charging  current  in  sine  wave  osc.  Vcc/(5*Ra) 

Isine  :  The  maximum  current  out  of  the  sine  wave  oscillator. 

Q  :  'Q'  for  filter  eq. 

:  'Q'  product  for  multiple  stage  filters. 

Ra,Rb  :  Charge  &  discharge  current  determining  resistors 

:  on  the  sine  wave  oscillator. 

Radj  :  Divider  resistor  for  sine  wave  osc  output. 

Rdiv  :  Thermistor  series  divider  resistor. 

Rfl  &  Rf2  :  Feedback  resistors  for  the  IP  filters  UAF41. 

Rfb  :  Feedback  resistor  of  coaiductance  Ist  anpllfrier. 

:  Determines  the  range  of  the  circuit. 

Rg  :  Gadn  resistor  for  the  IP  filter  UAF41. 

Rin  :  RMS  voltage  at  input  of  RMS — >DC  converter  in  conductance 

R1  :  Sine  wave  osc  load  resistor  or  Thermistor  divider  resistor 

Rout  :  DC  voltage  at  output  of  RMS — >DC  converter  in  conductance 

I^b  :  Probe  resistance  (  Conductance  ) . 

Rg  :  'Q'  determining  resistor  for  the  IP  filters. 

Rthm  :  Tenperature  probe  resistance. 

^>.gr.  :  Specific  gravity 

Teap  :  Tenperature  of  Oxvironnent . 

Vin  :  RMS  voltage  at  input  of  RMS — >DC  converter  in  conductance 

VGmax  :  Maximum  desired  sine  wave  voltage  at  input  to  Ist  anplifie 

VGtaiin  :  Minimum  desired  sine  wave  voltage  at  input  to  1st  anplifie 

VOrange  :  VCmin-VCmax 

Vout  :  DC  voltage  at  output  of  RMS — >DC  converter  in  conductance 

Vslg  :  RMS  voltage  measured  at  the  conductance  probe. 

Vsine  :  The  RMS  voltage  output  of  the  sine  wave  oscillator. 


tdSki. 


D«c  27, B4 
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INITIAL  PROGRAMMING  SETUP:  BAIN-1  OFFSET-O 
REMEMBER ! 

For  Btraight  binary  o-f-Taat  (0-10VDC>,  INVERT  MSBit 
From  -factory  addrass  decoding  -  BOH  —  Jumper  #1 


1.)  ZERO  TEMPERATURE  DAC 

A.  )  Zero  MSByte  1000  OOOO  1000  0100 

B.  )  Activate  MR.  (Pin  3) 

C. )  Zero  LSByte  OOOO  OOOO  lOOO  0101 

D.  )  Activate  WR.  (Pin  2) 

E. )  Xfer  data  to  DAC  output  XXXX  XXXX  lOOO  0110 

F.  )  Activate  WR.  (Pin  1) 


2. )  SET  COND.  PROG.  6AIN-1 

A. )  Set  Gain  of  Cond  Amp.  XXXX  XXOO  lOOO  OOll 

B. )  Activate  WR. 


3. )  ZERO  CONDUCTANCE  DAC  Db 

A.  )  Zero  MSByte 

B.  >  Activate  WR.  (Pin  3> 

C.  >  Zero  LSByte 

D.  )  Activate  WR.  (Pin  2) 

E.  )  Xf t»r  data  to  DM2  output 

F.  )  Activate  WR.  (Pin  1> 


7654  3210  Ad  7654  3210 

1000  OOOO  10<X>  (XK)0 

OOOO  (X)(X>  1(X>0  <X)01 

XXXX  XXXX  1(XX>  (X>10 


4. >  SET  TEMP  PROG.  GAIN-1 

A. )  Set  Gain  of  Temp.  Amp.  XXXX  XX(K>  lOOO  Olll 

B.  )  Activate  WR. 

EXAMPLE:  ASSUME  TEMP-1. 25  volte 

SET  OFFKT  TO  ZERO,  INCREASE  BAIN  TO  10 

Aeeummed  addrees  decoding  of  BOH 

1.)  Set  MSByte  OF  Temp  DAC  to  1.2SV 
LSByte  mill  remain  eame. 


Db  7654  3210  Ad  7654  3210 

A. )  Set  DAC  to  1.25V  1010  OOOO  lOOO  OlOO 

B.  )  Activate  WR. 

2.)  Set  Temp  Amp  Sain  to  10 

A. )  Set  Bain  XIO-  01  XXXX  XXOl  lOOO  Olll 

B.  )  Activate  WR. 
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FRCXSRAhMING  TABLES 


ZABI£  1 

ADDRESS  SELECT 


OCT  23.84 


Switch  A7  A6  AS  A4  A3  A2  A1  AO 

1  1000  OXXX  $BOtX 

2  1000  IXXX  $8BtX 

3  1001  OXXX  $90tX 

4  1001  IXXX  $98-t-X 

5  1010  OXXX  $A04-X 

6  1010  IXXX  $A8tX 

7  1011  OXXX  $BOtX 

8  1011  IXXX  $B8tX 

TABLE  2 

DAC/GAIN  SELECT 


A2  A1  AO 

Sel  Cond  DAC  MS  Byte  0  0  0 

Sel  Coni  DAC  LS  Byte  0  0  1 

Output  Cond  Dec  010 

Cond  Anpllfier  Gain  Oil 

Sel  Cond  DAC  MS  Byte  10  0 

Sel  Cond  DAC  LS  Byte  10  1 

Output  Cond  Dac  110 

Cond  Aiqpllfier  Gadn  111 

TABLE  3 

AMPLIFIER  GAIN 

Av  D1  DO 
10  0 
10  0  1 

100  1  0 

1000  1  1 


Note  1 
See  Table  3 


Note  1 
See  Table  3 


Note  1 

Data  viritten  to  this  register  is  ignored.  The  effect  o 
write  to  this  address  is  to  transfer  the  LSByte  and  MS 
a  16  bit  register  within  the  DAC.  A  change  in  the  outp 
only  occurs  after  a  write  to  this  register. 
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FIG  5A  BALANCED  DIVIDER  CIRCinT 

NOV  6,84  TABLE  OF  TERMISTC»%  RESISTANCES  VS  TEMPERATURE 
TDff*B^TURE  RANGE:  -5.0  — >  +25.0  Degrees  C 

Ref.  THERMQMETRICS  T.C.  Manual  181-D  page  A-2 

R25:  Zero  pcwer  res.  of  T.C.  at  cadibrated  tenp  25  deg.Ct273.15Kelvln 

RT:  Zero  poMer  res.  of  T.C.  at  tenp.  t(C)t273.15  Kelvin 

Rib:  Zero  power  resistance  at  tenp.  To;ToBt(C)-i-273.15 

T;  Absolute  teoperature  T:T*t(C)+273.15 

B:  Beta  of  T.C.  ...  Function  of  T.C  materiad. 

R25/R125:  Ratio  of  RT  @  2SDeg  to  RT  @  12SDeg 

RD1.FD2:  Divider  resistors  for  badanced  iiput  divider  network. 


FOR: 

10000  RTobR25=10K  Ohms 

10000  RDl  THM  divider  res . . .  positive  side 
10000  RD2  THM  divider  res . . .  negative  side 
3540  B^eta  is  Kelvin  tenp.  const,  for  mater iad  "A" 

19.8  R25/R125sl9.8 

273.15  Kelvin  ref. 

298.15  To*25+273.15 

2.5  "es;  Thermistor  divider  source  voltage 
RT  *  RTo  exp  {B(l/T  -  1/To)}  eo(T)  -eo;  eo(T)-es/[l+(RT/R) ] 

0.6  ■s;  Choosen  from  best  linearity 
curve  for  material  "A"  §  0  Deg  C 
15000  sRDl4RD2»R25/s ;  Resistance  needed  for  best 
0.002442  s  Data  System  Resolution  for  12b  Conv 


Rdl 

Rd2 

Rthm 

Degrees  C 

—6 

RT(C) 

K  Ohms 
39.55 

dRT 

K  Ohms 

View 

2.5 

Vhigh 

Voltage  Drop 

-5 

37.65 

1.903 

0.43367 

2.06633 

1.633 

-4 

35.85 

1.798 

0.44764 

2.05236 

1.605 

-3 

34.15 

1.700 

0.46169 

2.03831 

1.577 

-2 

32.54 

1.608 

0.47582 

2.02418 

1.548 

-1 

31.02 

1.521 

0.49001 

2.00999 

1.520 

0 

29.58 

1.440 

0.50424 

1.99576 

1.492 

1 

28.22 

1.363 

0.51849 

1.98151 

1.463 

2 

26.93 

1.291 

0.53276 

1.96724 

1.434 

3 

25.70 

1.223 

0.54702 

1.95298 

1.403 

4 

24.54 

1.159 

0.56126 

1.93874 

1.377 

5 

23.44 

1.099 

0.57546 

1.92454 

1.349 

6 

22.40 

1.043 

0.58961 

1.91039 

1.321 

7 

21.41 

0.989 

0.60370 

1.89630 

1.293 

8 

20.47 

0.939 

0.61771 

1.88229 

1.265 

9 

19.58 

0.892 

0.63163 

1.86837 

1.237 

10 

18.73 

0.847 

0.64544 

1.85456 

1.209 

11 

17.93 

0.805 

0.65913 

1.84087 

1.182 

12 

17.16 

0.765 

0.67270 

1.82730 

1.155 

13 

16.44 

0.727 

0.68613 

1.81387 

1.128 

14 

15.74 

0.692 

0.69941 

1.80059 

1.101 

15 

15.09 

0.658 

0.71252 

1.78748 

1.075 

16 

14.46 

0.626 

0.72547 

1.77453 

1.049 

17 

13.86 

0.596 

0.73825 

1.76175 

1.024 

18 

13.30 

0.568 

0.75084 

1.74916 

0.998 

19 

12.76 

0.541 

0.76324 

1.73676 

0.974 

20 

12.24 

0.515 

0.77544 

06 

1.72456 

0.949 

21 

11.75 

0.491 

0.78744 

1.71256 

0.' 

22 

11.28 

0.468 

0.79923 

1.70077 

0. 

23 

10.83 

0.447 

0.81082 

1.68918 

0. 

24 

10.41 

0.426 

0.82218 

1.67782 

0. 

25 

10.00 

0.407 

0.63333 

1.66667 

0. 

FIG  5A 

NOV  6 

,84 

DATA  SYSTEM 

RESOLUTION  V. 

s.  TEMP 

E  R  A  T  U  R  E 

@  OP  AMP  GAI 

N  S  : 

BASED  ON  lOVDC  F.S.  /  4095 

Degrees  C 

Gain  XI 

Gain  XIO 

Gain  XlOO 

Gain  XIOOO 

millldeg  C 

nC 

uC 

uC 

-5 

11.4376 

1 . 1438 

114.3756 

11.4376 

-4 

11.5111 

1.1511 

115.1109 

11.5111 

-3 

11.5716 

1 . 1572 

115.7163 

11.5716 

-2 

11.6191 

1.1619 

116.1910 

11.6191 

-1 

11.6535 

1.1653 

116.5348 

11.6535 

0 

11.6748 

1.1675 

116.7480 

11.6748 

1 

11.6832 

1.1683 

116.8316 

11.6832 

2 

11.6787 

1.1679 

116.7872 

11.6787 

3 

11.6617 

1.1662 

116.6167 

11.6617 

4 

11.6323 

1.1632 

116.3228 

11.6323 

5 

11.5908 

1.1591 

115.9083 

11.5908 

6 

11.5377 

1.1538 

115.3769 

11.5377 

7 

11.4732 

1.1473 

114.7322 

11.4732 

6 

11.3979 

1 . 1398 

113.9786 

11.3979 

9 

11.3121 

1.1312 

113.1206 

11.3121 

10 

11.2163 

1.1216 

112.1629 

11.2163 

11 

11.1111 

1.1111 

111.1107 

11.1111 

12 

10.9969 

1.0997 

109.9691 

10.9969 

13 

10.8744 

1.0674 

108.7436 

10.8744 

14 

10.7440 

1.0744 

107.4397 

10.7440 

15 

10.6063 

1.0606 

106.0630 

10.6063 

16 

10.4619 

1.0462 

104.6191 

10.4619 

17 

10.3114 

1.0311 

103.1137 

10.3114 

18 

10.1552 

1.0155 

101.5524 

10.1552 

19 

9.9941 

0.9994 

99.9407 

9.9941 

20 

9.8284 

0.9826 

98.2641 

9.8284 

21 

9.6588 

0.9659 

96.5679 

9.6588 

22 

9.4857 

0.9486 

94.8574 

9.4657 

23 

9.3098 

0.9310 

93.0976 

9.3096 

24 

9.1313 

0.9131 

91.3134 

9.1313 

25 
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OCT  31,84 

DESIGN  CALCDIATIQNS  FOR  BURR-BRCMN  11AF41  FILTER 
FUNCnCN  -  2  POLE  BUTIEIM3RTH  L.P. 


INVERTHIG 

2  >  #  of  poles. 

1  «  fn 
0.70711  ■  Q 

0.70711  <B  Qp  «  Q  for  fo*Q  products 
less  than  1E5  Hz 

25  >  fo(l) 

SO  -  fo(2) 

100  >  fo(3) 

B  «  Eqpiatlon  for  fo<8KHz 
Rfl«4lf2:  Feedback  Resistors 
Rfl-Rf2  -  5.033E7/(fo) 

Alp:  Low-Pass  Gain 
1  «  Alp 

Abp:  Banc^iass  Gain 
0.223768  «  Abp  «  (qp/3.16)«Alp 
Rg:  Gain  Resistor 
49928  -  Ret  -  (1.58E4*Qp)/Abp 
Rq:  Q  Resistor 

29682.19  >  Rq  >  5E4/[ (3.48*Qp)4-Abp-l] 

Redd:  Added  Res.  needed  for  Eq"B” 
5490  >  Radd  (Res  adkled:  pins  12&13} 

Cext:  Added  to  lower  Rf  values 
l.OOE-07  -  Cext  (.lufd) 

Rf  (new) :  New  Rfl  &  Rf2  be  to 

the  addition  of  Cext. 

Rf  ( new)  »^  ( old)  •  [  lE-9/ ( lE-94Cext )  ] 


fo  Rfl,Rf2  Cext  Rf(new)  Rg  Rq 

Kz  Ohns  Farads  Ohns  Ohms  Ohms 

25  2013200  1.00E-07  19932.67  IMV  49928  29682.19  Ohms 

50  1006600  1.00E-07  9966.336 

100  503300  l.OOE-07  4963.168 

S^ame  Rf  values 
INV  &  NON-INV 
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OCT  31,64 

CBSIGN  CALCULATIGNS  FCR  BURR-BRCMN  UAF41  FILTER 
FUNCnCN  -  2  POLE  BOTTERWDRIH  L.P. 

NDN-INVERTINS 

2  >  #  of  poles. 

1  -  fn 
0.70711  ■  Q 

0.70711  ■  Qp  ■  Q  for  fo*Q  products 
less  than  lES  Hz 

25  «  fo(l) 

50  -  fo(2) 

100  -  fo(3) 

B  ■  Equation  for  fo<6KRZ 
Rfl^^2:  Feedback  Resistors 
Rfl»Rf2  ■  5.033E7/(fo) 

Alp:  lOM-Paaa  Gain 
1  ■  Alp 

Abp:  Ban^)a88  Gain 
0.223768  -  Abp  «  (Q/3.16)«Alp 
Rg:  Gain  Resistor 
223444,7  -  Rg  »  (5E4*Q)/(Abp*Qp) 

Rq:  Q  Resistor 

40421.22  -  Rq  -  5E4/[3.48*Qp-(  (Abp*(i))/Q)-1] 

Redd:  Added  Res.  needed  for  Eq"B" 

5490  >  Radd  (Res  added:  pins  12&13) 

Cext:  Added  Ca^  to  ICMer  Rf  v^ues 
l.OOE-07  -  Cext  (.lufd) 

Rf  (new) :  New  Rfl  &  Rf2  due  to 

the  addition  of  Cext. 

Rf  (new)^(old) ‘[lE-a/dE-a+Cext)  ] 

to  Rfl,Rf2  Cext  Rf(nBw)  Rg  Rq 

Hz  Ohms  Farads  Otas  Ohns  OIns 

25  2013200  1.00E-07  19932.67  NW-IMV  223444.7  40421.22  ChoiS 
50  1006600  l.OOE-07  9966.336 

100  503300  l.OOE-07  4983.168 


Sane  Rf  values 
INV  &  NDN-INV 


TB<PERATURE  -  CONDUCTANCE 


Hav  12,84 


J1  Data  Txanslatlon  012745  20pin  Interccnnect 

J1 


CCmJCTANCE  SIG  CH  0  1 
TEMPERATURE  SIG  CH  1  3 
FFOH  {J2-P1)  CH  2  HI  5 
FROM  (J2-P3)  CH  3  HI  7 
(N/C)  Ch  4  HI  INPUT  9 
(N/C)  Ch  5  HI  INPUT  11 
(N/C)  Ch  6  HI  INPUT  13 
(N/C)  Ch  7  HI  INPUT  15 
(N/C)  ANA  CCM4DN  17 
FROM  (J2-P9)  EX  IRQ  19 


2  CCNDUCTANCE  SIG  GND 
4  TEMPERATURE  SIG  GND 
6  FRCM  (J2-P2)  CH  2  LO 
8  FROM  (J2-P4)  CH  3  LO 
10  Ch  4  LO  INPUT  (N/C) 
12  Ch  5  LO  INPUT  (N/C) 
41  Ch  6  LO  INPUT  (N/C) 
16  Ch  7  LO  INPUT  (N/C) 
18  AMP  CCMCN  (N/C) 

20  DIGITAL  GND  (N/C) 


J2  Sensor  I/O  lOpln  Connector 

J2 


Ch2  HE  INPUT  (J1-P7)  1 
Ch3  HI  INPUT  (J1-pT)  3 
Conciuctance  Pb  ->  5 

Thn  Pos.  Inpit  ->  7 

EXT  TRIG  (J1-P19)  9 


2  Ch2  LO  INPUT  (J1-P8) 
4  Ch3  LO  INPUT  (Jl-P^) 
6  <-  Ccs)diictance  GND 

8  <-  nim  Neg.  Input 

10  ->  ANALOG  GND 


O 


c* 
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riD-fll6S  258  8  HIGH-RESOLUTION  CLUSTER  OF  OCEHNOGRRPHIC  INSTRUNENTS  2/2 

FOR  BOUNOflRV  LflVER  HEASUREHENTS  UNDER  ICEIU)  FLON 
INDUSTRIES  INC  KENT  Ufl  H  LIU  ET  HL  NOV  85 
UNCLASSIFIED  TR-TSE/ll-SS  N88814-82-C-87fi4  F/G  S/IA  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATtONAL  BUREAU  OF  STANDARDS -1 963-A 


PTC- 1 

'KJCTANCE VOLIAGL  Novii.bf 


84 


OGMPQNENTS 


R1,R2 

lOK  1% 

MK132  C1,C4  .1  Monocap 

R3 

50K  1% 

R4 

30K  1% 

02,3,5,6,7,8,13,14,15,16 

R5 

100  1% 

.1  Poly 

R6 

200  Baums  pot  Model  3059P 

R7 

SK  1% 

05,9,10,17,18 

R8 

3QK  1% 

.1  Mono 

R9 

50K  151$ 

RlO.Rll 

lOK  1% 

C11,C12  .01  Poly 

R12«R13 

lOK  Carbon  55|$  C19,C20  68pfd  Mica 

R14 

82K  151$ 

R15,R16 

lOK  151$ 

021,22,23,24 

RNl.2,3. 

4  IQK  151$ 

10  i  35V 

694-3-RlOK-F  8PIN  DIP 

Q1 

AD580 

Vref.  2.5VDC  C25  4700pf 

U1 

74LS02 

NOR  GATE,  LATCH  DB  WHEN  ADDR  &  Ml  ARE  TRUE 

U2 

74LS221 

HR  PULSE  STRETCHER  MC»IO  FOR  FGA*S 

U3,U4 

74LS138 

ADDRESS  DECODING 

US 

74LS373 

DB  LATCH 

U6 

ICL8038 

PREC.  SINE  HAVE  GENERATOR 

U7,U8 

DAC709BG 

16  BIT  D/A,  OFFSETS 

D9,U10 

AD536 

RMS  — >  DC  CONVERTERS 

Ull 

AD534 

DIVIDOl,  Isig/VrefsG  (conductance) 

012,13,20,21,22 

OPAlll  PREC  CP  AMP 

U14,U15 

PGA200BG 

PROG  GAIN  AMPS.  X1,X10,X100,X1000 

U16,U17 

If411CN 

MED  SPEED  OP  AMPS 

U18,U19 

UAF41 

L.P.  FILTERS 

J1 

20  PIN  HEADER  CGNNECTC3R,  MALE.  TO  DT2745  A/D  ED 

J2 

10  PIN  HEADER  CONNECTOR,  MALE.  INPUTS 

WV  12,84 


PTC-1  CIRCUIT  BD 


5  :  Minimum  #  of  boards  needed. 
450  :  Total  #  of  holes. 

« 0.003  :  Hole  tolerance,  all  holes. 
29.25  :  Board  area  (  sq.  in.  ) 

0.062  :  Board  thidmess. 


SPECIAL  INSTRUCTIONS: 

BRIGHT  ACID  TIN  TRACES. 
GOID  PLATED  FINGERS. 
SOIDER  MASK. 

LEGEND. 


Pod  Master 


